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Abstract 

The nitrogen-vacancy (NV) colour centre in diamond is an im- 
portant physical system for emergent quantum technologies, includ- 
ing quantum metrology, information processing and communications, 
as well as for various nanotechnologies, such as biological and sub- 
diffraction limit imaging, and for tests of entanglement in quantum 
mechanics. Given this array of existing and potential applications 
and the almost 50 years of NV research, one would expect that the 
physics of the centre is well understood, however, the study of the NV 
centre has proved challenging, with many early assertions now be- 
lieved false and many remaining issues yet to be resolved. This review 
represents the first time that the key empirical and ab initio results 
have been extracted from the extensive NV literature and assembled 
into one consistent picture of the current understanding of the centre. 
As a result, the key unresolved issues concerning the NV centre are 
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identified and the possible avenues for their resolution are examined. 
Keywords: diamond, color center, nitrogen-vacancy 

1 Introduction 

The detection of single negatively charged nitrogen- vacancy (NV~) colour 
centres in 1997 [1] marks a critical point in the evolution of diamond based 
quantum technologies. Although observations of ensembles of NV~ centres 
were routinely performed prior to 1997, the detection of single centres soon 
enabled demonstrations of photostable single photon generation [21 [31 H], 
which highlighted the NV~ centre for implementation in quantum optical 
networks, as well as demonstrations of optical preparation and readout of 
the centre's electronic spin [5l [6] , which identified the NV~ centre as a pos- 
sible solid state spin qubit suitable for quantum information processing and 
quantum metrology devices. Following these demonstrations, the growth of 
research into the NV~ centre and the development of applications of the cen- 
tre have been incredibly rapid and a number of important milestones have 
been reached. Although the NV centre is also found in the neutral charge 
state (NV°), NV° has not yet been employed in any demonstrations akin to 
NV~ and there exists only very speculative arguments [7] concerning its suit- 
ability as a spin qubit (due to no reported observations of optical readout). 
Consequently, the motivations for studying NV° have thus far been restricted 
to obtaining insight into NV~. 

The NV~ centre has been employed in milestone room temperature demon- 
strations of quantum registers built upon the NV~ electronic spin and prox- 
imal N and ^^C nuclear spins [SI [H]. At low temperature (< 10 K), spin- 
photon entanglement between the ground state spin of a single NV~ cen- 
tre and the polarisation of an emitted photon has also been demonstrated 
|10j . In terms of scalability, demonstrations of NV~-NV~ spin coupling 
[TT] have been performed and important steps towards photonic coupling 
have been made [12]. Beyond quantum information processing, there ex- 
ist important applications in room-temperature nanoscale magnetometry 
[laiHIiailSlIIZllIHllIHlEOlEIlEaEn]^ blo-magnetometry [21] , electrometry 
[25] and decoherence microscopy [261 EZl ESj using the NV~ centre. Adding 
to the already impressive list of metrology applications, the NV~ centre has 
also been recently proposed as a sensitive nanoscale thermometer [211 ED] • 

Since the demonstrations of the N V~ centre as a solid-state spin qubit and 
a nanoscale multi-sensor, the centre has been used to perform experiments ad- 
dressing fundamental issues in quantum and condensed matter theory. These 
experiments include Bell state violation [31] , magnetic resonance beyond the 
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rotating frame approximation [32], detection of the Meissner effect [33] and 
quantum process tomography [31|. Recently, the NV~ centre has been pro- 
posed as an ideal system to study mesoscopic quantum entanglement between 
an NV~ spin qubit and the zero-point vibrations of a magnetic micromechan- 
ical cantilever, or alternatively, the use of a cantilever to entangle NV~ spin 
qubits separated by mesoscopic distances [33, 13S] • As a photostable single 
photon source, the NV~ centre has been utilised to demonstrate quantum 
cryptography [371 [38] and the generation of single surface plasmons on silver 
nanowires [39] . The ability to tune the optical transition of the NV~ centre 
using strain and electric fields has also led to the centre being proposed for 
tests of cavity quantum electrodynamics [121 HOl IH] and photonic devices 
such as Q-switching photonic modules [121 113] • 

In the almost 5 decades of NV research, many properties of both charge 
states have been established. However, the understanding of these properties 
has not been without contention and significant revision and there remain a 
number of properties that are not well understood. The notable past issues 
concerning the NV~ centre have been the identification of its charge state, 
the spin multiplicity of the ground electronic state as well as the intermedi- 
ate electronic states involved in the non-radiative decay from the optically 
excited electronic state to the ground electronic state that enables the op- 
tical preparation and readout of the centre's electronic spin. The ongoing 
issues concerning the NV~ centre include the detail of the aforementioned 
non-radiative decay, anomalous aspects of the centre's vibronic structure, 
the Jahn- Teller effect in the centre's optical transition and temperature vari- 
ations of the centre's properties. Other important aspects of the NV centre 
that are not well understood are the factors that govern the relative charge 
state concentrations, the mechanism of photoconversion between the charge 
states and the associated spectral diffusion of the NV~ optical transition. The 
theoretical methods that have been employed to understand the NV centre's 
properties have been the semi-empirical molecular model of deep level defects 
in semiconductors [lOlliailSlllHlllZlllHlllSlEOlEIlEaE^ tmtio cal- 

culations [71|5ai55l[Ml[571|58l|59l[60l[6ll|621l6^ 

These theoretical methods have complementary strengths and weaknesses, 
which necessitate their combined application to address the ongoing issues. 
Clearly, the improved understanding gained from the resolution of these on- 
going issues will enable enhanced implementation of the NV centre in its 
many impressive applications. 

Both the rapid expansion of the field of NV research and the need to 
resolve the ongoing issues motivate this detailed review of the properties of 
the NV centre in which we identify and define ongoing issues and discuss 
avenues to their resolution. In addition we aim to provide a comprehensive 
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and self-consistent reference for those both new and old to NV research, which 
is free of any past interpretations that are now believed false. In the following 
subsections, the properties of the NV centre will be firstly introduced in order 
to orientate those new to the centre and to also facilitate the subsequent brief 
introduction to the quantum technology applications of the centre that have 
motivated its intense investigation. The final subsection presents the scope of 
this review, where the domain of the NV literature considered in this review 
is defined and the method of the review is discussed. 

1.1 Introduction to the NV centre 

The NV centre is a point defect in diamond with C^^ symmetry consisting 
of a substitutional nitrogen - lattice vacancy pair orientated along the [111] 
crystalline direction [71] (refer to figured]). The centre may be found as an 'in 
grown' product of the chemical vapour deposition (CVD) diamond synthesis 
process [72] or as a product of radiation damage and annealing [73] or ion 
implantation and annealing [71] in bulk and nanocrystalline diamond. The 
centre is known to exist in negative (NV~) [HI and neutral (NV'') [75] 
charge states. The identifying features of NV~ and NV" are their optical 
zero phonon lines (ZPLs) at 1.945 eV (637 nm) [73] and 2.156 eV (575 nm) 
[76], respectively, and associated vibronic bands that extend from their ZPLs 
to higher/ lower energy in absorption/emission (see figured] for examples of 
the emission vibronic bands). An additional infrared ZPL at 1.190 eV (1042 
nm) that is only observable during optical illumination (with energy > 1.945 
eV) has been associated with the NV~ centre [77]. The sharp ZPLs and 
well defined vibronic bands of NV° and NV~ indicate that the optical transi- 
tions occur between discrete defect levels that are deep within the diamond 
bandgap, such that the continua of valence or conduction band levels are 
not involved in the optical transitions. The NV centre is thus known as a 
deep- level defect in diamond [lU HS] . A schematic of the known aspects of 
the electronic structures of NV~ and NV° are depicted in figure [21 

The NV~ centre may also be identified by a zero field magnetic reso- 
nance at ~2.88 GHz [H| B5] (see figure Ej). This magnetic resonance occurs 
between the rris = and = ±1 spin sub-levels of the spin triplet ground 
state and can be detected by either conventional electron paramagnetic 
resonance (EPR) [HI HH] or optically detected magnetic resonance (ODMR) 
[78] techniques. A further ODMR at ~1.42 GHz [791 [80] is observed at room 
temperature (see figure [3]) and is attributed to a zero field splitting of the 
spin triplet excited state that is analogous to that of ^^2. The magnetic 
resonances of and behave as would be expected of triplet spins in 
a trigonal crystal field and are characterised by the approximately isotropic 
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(a) (b) 

Figure 1: (a) Schematic of the nitrogen- vacancy centre and diamond lattice 
depicting the vacancy (transparent), the nearest neighbour carbon atoms 
to the vacancy (green), the substitutional nitrogen atom (brown), and the 
next-to- nearest carbon neighbours to the vacancy (black), (b) Normalised 
emission spectra of an ensemble of NV~ and NV° centres at 10 K for different 
excitation (532 nm) powers as indicated [83]. The photoconversion between 
the charge states is apparent in the increase of the emission of NV° relative 
to NV~ with increasing excitation power. 



electron g-factors ~2.0028 [SI HS] and ~ 2.01 [79] respectively. The mag- 
netic resonances of [SI |82] and [791 EO] also exhibit weak interactions 
with strain and electric fields. 

At low temperatures (<10 K), the excitation spectrum of the NV~ opti- 
cal ZPL (refer to figure [3]) reveals a fine structure of the ZPL that is highly 
dependent on strain [H] and electric fields [85] . The fine structure of the ZPL 
implies a fine structure of that is much more complicated than that sug- 
gested by the associated ODMR at room temperature. The low temperature 
fine structure of '^E and its observed dependence on strain [SI] is depicted 
in figure m (analogous results have been obtained using electric fields [53]). 
It is clear that the presence of strain and/or an electric field splits the ^E 
fine structure into two branches; the lower branch undergoes two level anti- 
crossings (as marked in figure H]) that mix the lower fine structure states to a 
degree dependent on the strain/electric field, whereas the upper branch does 
not experience any level anti-crossings and the upper fine structure states 
remain approximately unmixed [85] . The mixing of the lower fine structure 
states results in optical spin-flip transitions when the lower branch is excited. 
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Figure 2: Schematics of (a) known aspects of the NV~ electronic structure 
|l6] and (b) known aspects of the NV° electronic structure [86]. The observed 
optical [731 ESj and infrared [77] zero phonon lines (ZPLs) are as indicated. 
Note that the relative energies of the triplet (^E, ^^2) and singlet states {^Ai, 
^E) of NV" and the doublet (^A, ^E) and quartet states (^^2) of NV° are 
currently unknown. 



whilst the unmixed upper fine structure states maintain optical spin selection 
rules (AS* = 0, Am^ = 0) and result in optical spin-conserving transitions 
when the upper branch is excited [HSj. The collapse of the low temperature 
^E fine structure into a single broad ODMR at room temperature is thought 
to be due to phonon mediated orbital averaging over the ^E fine structure 
levels [87]. This orbital averaging process is depicted in figure |4] and has the 
consequence that individual centres with different ^E fine structures at low 
temperature (due to differing crystal strains) will have near identical '^E fine 
structures at room temperature. In addition to their role in orbital averaging, 
phonons are also believed responsible for the temperature dependent depo- 
larisation and broadening of the optical ZPL (see figure [5]) [88]. A thorough 
understanding of the phonon processes in ^E is yet to be obtained. 

The optical detection of the magnetic resonances of the NV~ centre is 
enabled by the differing fluorescence of the = and ±1 spin projec- 
tions [78]. As demonstrated in flgure El the spin dependence of the NV~ 
fluorescence can be explained by different non-radiative decay pathways for 
the ms = and ±1 spin projections from ^E to via intermediate dark 



6 



3 
CO 



c 
<u 
u 

(i> 
c 

E 

_3 

O 

■»-' 

o 




excited state 



1.42 GHz 



ground state 



:±1 



2.88 GHz 



1000 1500 2000 2500 

Frequency (MHz) 

(a) 



3000 




-10-5 5 

Laser frequency (GHz) 

(b) 



Figure 3: (a) Room temperature ODMR spectra [SU] of a single NV centre 
exhibiting the zero field splittings of the = and rris = ±1 spin sub- levels 
of (~ 2.88 GHz [mis]) and (~1.42 GHz [79]). (b) Low temperature 
(4 K) excitation spectra [H] of different single centres exhibiting the highly 
strain dependent NV~ optical ZPL fine structure. 



states (believed to be Mi and ^E) gZl [891 190]. As these non-radiative de- 
cay pathways compete with fiuorescent decay, the spin projection with the 
weakest non-radiative decay out of ^E will appear brightest and vice versa. 
The observed optical enhancement of the EPR signal strength implies 
that the probability of finding the centre in a particular spin projection is 
polarised by optical excitation, a process known as optical spin-polarisation 
[HI H5] . The combination of a number of observations has led to the con- 
sensus that the = spin projection is the bright projection and also the 
projection into which the centre is optically polarised [HH [221 [H31 [M]- It is 
currently believed that the non-radiative decay from the ^E spin sub-levels 
to the dark states is slower for = than for = ±1 and that the 
non-radiative decay out of the dark states to the spin sub-levels is sim- 
ilar for both spin sub-levels, thereby leading to a net polarisation into the 
nis = sub- level [TUl [13 [HSl [HSl [HZ] • Through observing the time-dependence 
of the centre's fiuorescence (refer to figure [S]), it is clear that the process of 
optical spin-polarisation takes a few optical cycles; the initial fiuorescence 
is representative of the initial spin state of the centre and spin-polarisation 
is achieved once the optical steady state has been reached. Consequently, 
the optical dynamics of the NV~ centre can be used to optically readout the 
centre's spin state and to simultaneously spin-polarise the centre into the 
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Figure 4: (a) The strain dependence of the fine structure observed at 
low temperature (4 K) [SI]. The level anti-crossings of the lower branch 
(Ey) are circled, (b) Schematic of the observed low temperature ^E fine 
structure at low (left) and high (right) strain together with the observed 
room temperature fine structure (centre) [17]. The spin-projection (m^ = 0, 
±1) of each fine structure level is as indicated. The phonon transitions that 
are proposed to give rise to orbital averaging are indicated by arrows |46j . 



nis = spin projection Specifically, optical readout is conducted by 

comparing the integrated fiuorescence of the spin state to be read out with 
calibration measurements of the integrated fiuorescence corresponding to the 
spin being prepared in the = and ±1 spin states [98] . 

Differing from NV~, NV° does not have detectable magnetic resonances 
associated with its spin doublet ground (^E) and excited (^A) states. The 
absence of a NV° ground state magnetic resonance is believed to be due to the 
presence of the dynamic Jahn- Teller effect [86]. It was recently discovered 
that NV° does, however, exhibit an EPR signal under continuous optical 
illumination [SB]. This EPR signal has been attributed to a long lived spin 
quartet excited state i^A2) that is populated by non-radiative decay from the 
optically excited "^A [86]. No ODMR or optical readout of the NV*^ quartet 
spin state have been demonstrated to date. Yet, the EPR signal strength 
has been observed to increase with optical excitation intensity, indicating an 
optically dependent spin-polarisation of ^ [SS] , albeit via a very different 
process to that which occurs in NV^. 
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Figure 5: (a) Temperature dependence of the fluorescence polarisation visi- 
bility {Ix — ly) / {Ix + W) for different single NV~ centres (assorted coloured 
dots) and different excitation polarisations (pump X/Y as indicated) [88] . 
The departure from 100% visibility signifies phonon transitions between the 
fine structure states that act to depolarise the fluorescence [88]. (b) 
A series of resonant excitation spectra scans of a single NV~ centre at 10 K 
demonstrating the spectral diffusion of the optical ZPL (dark spots) in terms 
of detuning from 1.945 eV over time |88j . 



Regardless of the competing non-radiative processes in the optical cycles 
of NV~ and NV*^, both charge states have sufficiently strong fluorescence to 
be detected as single centres using scanning confocal microscopy [T] (see figure 
[7] for an example). Indeed, both centres have been identified as single photon 
sources [21 [3l IH [99], as demonstrated by the observed photon statistics in 
figure [3 The detection of single centres implies certain characteristics of the 
NV~ and NV° optical cycles: large optical absorption cross-sections, short 
excited state lifetimes, high quantum yields, and no efficient shelving in a long 
lived dark state [T]. Furthermore, both charge states have been observed to be 
extremely photostable under off-resonance excitation (typically 2.32 eV (532 
nm)) [H [31 m [99] , with no evidence of photobleaching except in extreme condi- 
tions [lOOj . However, photoconversion between the charge states and photo- 
blinking under resonant optical excitation has been regularly observed (refer 
to figure [T|) [21 1101] , but it is not currently understood. The photoconversion 
efficiency has been observed to be dependent on optical excitation energy and 
intensity [IZl [831 [971 (M [M [M [I06l [M] as well as material proper- 
ties and the method of NV formation [JOH [I09l [IIOl nH [112]. The spectral 
diffusion of the NV~optical ZPL (refer to figure [S]) is also believed to be re- 
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(a) (b) 

Figure 6: (a) Schematic of the NV~ electronic structure including the 
and (room temperature) fine structure [l6]. The optical and infrared 
transitions are denoted by solid arrows and the weak (— ■ — ) and strong (- 
- -) non-radiative transitions to, between and from the intermediate singlet 
states are denoted by dashed arrows, (b) Examples of the time-dependent 
fiuorescence Im^ for pure initial spin states {rris = 0,±1) [98]. The data 
was collected from a single centre using off-resonant optical excitation and 
accumulated over many measurement repetitions. The contrast in the flu- 
orescence of the pure spin states is demonstrated in the lower half of the 
plot. 



lated to photoconversion and the ionisation of electron donors in the vicinity 
of the NV" centre during optical excitation [881 EH IMl [HB M EM EH- 
The focus on individual NV~ centres stems from the large inhomogeneous 
strain broadening of the optical ZPL and signiflcant spin cross-relaxation 
effects present in high density ensembles. These ensemble effects prevent 
the resolution of the NV~ low temperature optical ZPL flne structure and 
broaden ODMR lines, thereby reducing the visibility of many of the centre's 
properties. The evolution towards the study of individual NV centres has 
been enabled by advances in single molecule spectroscopy [115^ 1116] , NV 
formation techniques and diamond synthesis. The advent of CVD synthesis 
techniques provided the means to control the size and quality of the diamond 
crystal [13]. The employment of nitrogen ion implantation in combination 
with an annealing step has facilitated the (non-deterministic) creation of 
single NV centres with nanoscale accuracy |117] . Apart from bulk diamond. 
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Figure 7: (a) Example room temperature confocal image [118] of single NV~ 
centres, (b) Example photon autocorrelation function |118j obtained from a 
single NV~ centre at room temperature. The antibunching dip below 0.5 at 
r = is indicative of a single photon source and confirms the presence of a 
single NV~ centre |118j . 



single NV centres have been fabricated in nanodiamonds as small as ~5 
nm in diameter |100] . Consequently, there exists a diverse range of crystal 
environments in which the NV centre can be formed and their effects on 
the centre investigated. Indeed, the advances in spectroscopy and material 
science have provided a suite of tools to study the NV centre and develop its 
apphcations. 

1.2 Overview of the NV in the quantum technologies 
context 

The applications of the NV~ centre each utilise one or a combination of the 
following five key properties of the centre: 

1. A bright photostable optical transition that is suitable for the detection 
of individual centres and single photon generation 

2. An optical ZPL fine structure that is dependent on electric, magnetic 
and strain fields at low temperature, but is approximately electric and 
strain field independent at room temperature 

3. A magnetically resonant and controllable ground state electronic spin 
that exhibits long coherence times and coupling to proximal electronic 
and nuclear spins 
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4. Optical spin-polarisation and readout of the ground state spin 

5. Flexibility and robustness in fabrication 

Some aspects of the above properties have already been discussed. Expand- 
ing on the second property in the context of the centre's apphcations, the 
strain and electric field dependence of the low temperature optical ZPL fine 
structure enables the transition energy of the ZPL to be tuned and spin- 
flip/conserving transitions to be selectively excited. The approximate strain 
and electric field independence of the room temperature fine structure al- 
ternatively offers some degree of uniformity in the fine structure between 
centres and simplifies the microwave control of the excited state spin of 
the centre. Whilst the magnetic resonance of the ground state spin has 
been discussed already, its coherence has not been addressed. Indeed, the 
ground state spin has the longest room temperature single spin coherence 
time (T2) of any electronic spin in a solid, being greater than 1.8 ms in 
highly engineered samples [13]. The long coherence time of the ground 
state spin enables its couplings with proximal electronic and nuclear spins 
in the lattice to be resolved and manipulated (see figure M for an example) 
El El [m [IISl EM IMl [Eg. The final property refers to the advances in 
material science that have developed a wide range of means to deliberately 
fabricate both individual centres and ensembles of NV~ centres in bulk and 
nanodiamond (refer to figure [S] for an example), as well as the observed ro- 
bustness of the first four properties to the fabrication techniques employed 

Each of the room temperature implementations of the NV~ centre as a 
spin qubit employ the same prepare-manipulate/interact- readout mode of 
operation (see figure (H]). An initial off- resonance optical pulse is used to pre- 
pare the ground state spin via optical spin-polarisation. This is followed by 
the application of static and microwave fields to manipulate the spin (and 
potentially radio frequency fields to manipulate coupled nuclear spins) and 
a period to allow the NV to interact with other spins and evolve. Finally, an 
off-resonance optical pulse is used to readout the ground state spin by mea- 
suring the integrated Stokes shifted fluorescence. It is clear that this mode 
of operation relies upon the performance of the optical spin-polarisation and 
readout mechanisms, which are currently only crudely understood. Indeed, 
the degree of ground state optical spin-polarisation is not consistently re- 
ported in the literature, with many different values ranging from 42-96% 
reported El [IDl IHl IMl EHl IS31 [Ml EM EM [I3D]. Since the ground state 
spin-polarisation represents the preparation fidelity of the qubit, it is vital 
to the operation of the qubit that it is well characterised. However, to date 
there has been no systematic study of spin-polarisation and its variation with 
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(a) (b) (c) 

Figure 8: (a) Schematic of the proximal nuclear spins (arrows) to the NV 
centre: the native nitrogen nuclear spin (^^N or ^^N) and the randomly dis- 
tributed ^^C nuclear spins that are naturally 1.1% abundant in diamond, (b) 
Example ODMR spectra [9] demonstrating the hyperfine structure due to 
interactions with two ^^C nuclear spins (no coupling to the N nuclear spin 
observable). Below is the corresponding hyperfine structure with the logical 
bits (0,1) representing nuclear spin projections (m/ = — +i) [H]. Electronic 
spin transitions observable in the ODMR spectra are denoted by solid arrows 
and the nuclear spin transitions that can be used for control are denoted by 
dashed arrows [9]. (c) Example spatial photoluminescence image of a co- 
planar waveguide (CPW) and regular structure of individual NV~ centres 
(white dots) created by N"*" ion implantation |117] . 



applied electromagnetic fields, strain, and temperature. A second significant 
problem is the readout contrast, which due to low collection efficiencies of 
current apparatus, is limited to a difference of a fraction of a photon be- 
tween the rris = and ±1 spin-projections |98 | 1131] . Consequently, a single 
qubit operation must be performed many times before sufficient shot-noise 
statistics are acquired to distinguish the spin-projections pS j 1131] . Recently, 
there have been reports of readout enhancement techniques using hyperfine 
coupling with N [98] and ^^C |131] nuclear spins, however it is quite possible 
that a fully developed understanding of the optical readout mechanism will 
provide alternate means to enhance readout more significantly. 

The use of the NV as a nanoscale sensor has been proposed for DC, AC 
and fiuctuating magnetic fields [HI [151 ESI ETJ 1132] . Each of the sensing ap- 
plications utilise one or more of three possible architectures (refer to figure 
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Figure 9: Schematic of the prepare-manipulate/interact-readout mode of 
operation of the NV~ spin qubit (hght blue circle with arrow). An ini- 
tial off-resonance (green) optical pulse is used to prepare the spin via opti- 
cal spin-polarisation. Control fields such as microwaves and static electric, 
magnetic and strain fields are applied to manipulate the spin and the other 
spins (other circles with arrows) that it is interacting with. Finally, an off- 
resonance (green) optical pulse is used to readout the spin via the integrated 
red shifted fluorescence. The readout pulse also prepares the spin for the 
next manipulate/interact step. 



[TOj) : a single NV~ centre in a nanodiamond attached to a scanning atomic 
force microscopy (AFM) cantilever, an array of individual NV~ centres in a 
bulk crystal substrate or in separate nanodiamonds, and an ensemble of NV~ 
centres in a bulk crystal substrate. The spatial resolution and sensitivity of 
these architectures are predominately governed by the proximity of the NV~ 
centre to the sample (as determined by the radius of nanodiamond or depth 
of the NV~ centre from the surface of the bulk crystal), the coherence time 
of the NV~ ground state spin and the number of NV~ centres in the sensing 
volume [T^[T5] . Consequently, given the inclusion of NV~ centres in nanodi- 
amonds, the ability to fabricate high density ensembles or arrays of centres 
with nanoscale tolerances in bulk diamond, and sub-diffraction limit imaging 
and magnetic field gradient techniques, the NV~ centre offers nanoscale spa- 
tial resolution in each of the three architectures [TU [15]. Isotopically pure 
diamond [IH] and dynamic decoupling techniques [221 11331 1134^ 1135] both 
yield effective spin coherence times of T2 > 1.8 ms and magnetic field sen- 
sitivities of single centres as high as ~4 uT/a/Hz [13]. Each of the three 
architectures employ the same mode of operation as the implementations of 
the NV~ centre as a spin qubit and are thus subject to the same polarisation 
and readout problems. However, additional problems arise in the centre's 
sensing applications because of the ground state spin's sensitivity to not 
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only magnetic fields, but to also electric and strain fields [251 1136] . Being a 
multi-field sensor, the NV~ centre must be optimised for each sensing mode, 
otherwise the sensitivity to one field will be inhibited by unwanted sensitivity 
to the other two fields |25l 1136] . A further added complexity is the temper- 
ature dependence of the zero field magnetic resonance of the ground state 
spin [22], which is currently not well understood and must be a consideration 
in the optimisation of the sensor. Due to the potential for significant tem- 
perature fluctuations when operating at room temperature, the temperature 
dependence of the NV~ spin will particularly limit the field sensitivity of the 
spin when operating at room temperature if it is not correctly optimised [29]. 
Alternatively, the NV~ spin may be optimised to sense the fluctuations in 
temperature and acts as a nanoscale high sensitivity thermometer [2^1 15U] . 



AFM cantilever 




(a) (b) 

Figure 10: (a) Schematic of one of the NV~ field sensing architecture employ- 
ing an NV~ centre in a nanodiamond attached to an AFM cantilever to scan 
a sample. As for the centre's qubit implementations, the field sensing im- 
plementations use optical preparation and readout and static and microwave 
fields for control of the NV~ ground state spin, (b) Example vector magnetic 
field measurement of the magnetic field induced by parallel current carrying 
wires using an ensemble of N V~ centres in a diamond substrate [20] . 



1.3 Scope of this review 

This review is the first time that the key empirical and ab initio results 
have been extracted from the extensive NV literature and assembled into 
one consistent picture of the current understanding of the centre. Given this 
objective, the scope of this review has been restricted to the intrinsic proper- 
ties of the NV centre and rather than the centre's various applications. There 
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are a number of reviews in existence that deal specifically with the centre's 
apphcations 1118^ 1137^ 1138] . Motivated by the numerous applica- 

tions of the negatively charged system NV~, it will be the primary subject of 
this review. However, as it is clear that NV° offers valuable insights into the 
properties of the NV~, NV° will not be ignored. The rapid advance of NV 
literature also necessitates the restriction of this review to works pubhshed 
prior to 2012. 

The structure of this review is as follows. The first section examines the 
identification of the physical structure of the NV centre and its various charge 
states by the initial NV studies. The following electronic structure section 
details the key studies that have provided the current electronic structures of 
NV~ and NV°, including the variations of the structures due to electromag- 
netic and strain fields. This section also discusses the key ab initio calcula- 
tions of the electronic structures. The vibronic structure section discusses the 
optical and infrared vibronic bands of NV~ , including the asymmetry of the 
optical bands, the recent work on the Jahn- Teller effect within the and 

electronic levels and the temperature variations of the centre's properties. 
The final optical and spin dynamics section discusses the extrinsic (i.e. the 
photoconversion mechanisms) and intrinsic (i.e. the spin-polarisation mecha- 
nism) aspects of the NV~ optical dynamics as well as the intrinsic coherence 
properties of the NV~ ground state spin. After the comprehensive review 
offered in the preceding sections, the ongoing issues in the understanding of 
the NV centre are identified and the possible avenues for their resolution are 
examined in the conclusion. 

2 Defect structure and charge states 

In 1965, du Preez observed the formation of the NV^ optical band in type 
lb diamond following radiation damage and annealing [73]. Remarkably, 
he speculated that the band was most likely generated by a substitutional 
nitrogen- vacancy pair, since type lb diamond contains a significant concentra- 
tion of single substitutional nitrogen (N^) defects and the process of radiation 
damage and annealing creates mobile lattice vacancies. Du Preez supported 
his argument by an observation of the correlated decay of the intensity of 
the optical band of the lattice vacancy (GRl) centre with the growth of the 
intensity of the NV~ optical band. The trigonal symmetry of the centre and 
the A-E electronic nature of the optical transition were established by the 
photoluminescence polarisation study conducted by Clark and Norris |139j 
and the subsequent uniaxial stress study conducted by Davies and Hamer 
|71j . The establishment of the centre's symmetry supported du Preez's spec- 
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ulation of a nitrogen- vacancy pair formed along the [111] direction. Davies 
and Hamer took du Preez's speculation a step further and concluded that 
this must indeed be the structure of the centre, using their analysis of the 
optical absorption and luminescence bandshapes. 

Davies and Hamer noted a slight asymmetry in the two bandshapes, which 
is most pronounced in a small splitting of the one-phonon peak in the absorp- 
tion band that does not occur in the luminescence band (refer to figure [TT] for 
a comparison of the absorption and luminescence bands). Such a vibronic 
splitting is typically the signature of the Jahn- Teller effect |140] . however, 
Davies and Hamer observed that after eliminating random lattice strains in 
the NV~ ensemble, the remaining depolarisation of the one-phonon peak of 
the luminescence band (typically taken to be secondary evidence of the Jahn- 
Teller effect |14U] ) was exceedingly small and they concluded that there was 
no significant Jahn- Teller effect. Instead, Davies and Hamer adopted an ex- 
planation of the vibronic band asymmetry based upon the elegant notion that 
the nitrogen nucleus tunnels between its substitutional site and the vacancy 
site. The tunneling model is analogous to, and no doubt drew inspiration 
from, the inversion tunneling of nitrogen in ammonia. Whilst the tunneling 
model was attractive in its elegance, it appears to have led to a fortuitous 
identification of the NV centre, given the recent empirical [881 1141[ I142j and 
ab initio [SU [601 ES EO] evidence in support of the Jahn- Teller explanation. 
The failure of the tunneling model and the readoption of the Jahn- Teller 
explanation will be discussed later in this review. 

Confirmation of the nitrogen-vacancy structure and the proposal of the 
negative charge state were provided by Loubser and van Wyk [HlllS] through 
their observation and interpretation of an EPR signal correlated with the 
NV~ optical band. The EPR spectra corresponded to that of a triplet 
state, whose presence Loubser and van Wyk explained by six electrons occu- 
pying two Ai (ai, a[) and degenerate E {cx, e^) molecular orbitals (MOs), 
formed from the dangling bonds of the nitrogen and carbon atoms surround- 
ing the vacancy, such that the two unpaired electron spins of the triplet state 
occupied the E MOs (see figure [12] for Loubser's and van Wyk's model). 
This most simple version of the molecular model was typical of the models 
used to explain EPR spectra of defects in semiconductors at the time [143] . 
The adoption of a molecular model of the NV centre has been supported by 
both EPR [861 HSOl MM and ab initio studies of NV" and NV°, which have 
confirmed that the NV MOs (ai, e^, Cy) within the diamond band gap are 
highly localised (> 72%) to the nearest neighbours of the vacancy. Note that 
the majority of ab initio studies place the lowest energy MO predicted by 
Loubser and van Wyk within the diamond valence band. 

In a chemical bonding picture, five of the six electrons were attributed 
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Figure 11: Comparison of the absorption (solid line) and luminescence (dot- 
ted line) bands of NV~ [7T]. Note that the absorption band has been reflected 
in energy about the ZPL at 1.945 eV in order to directly compare with the 
luminescence band. The ZPL present at 1.819 eV is due to another defect 
and not the NV centre |71j . 



to the electrons occupying the dangling bonds of the nitrogen and carbon 
atoms surrounding the vacancy [HlllS]. The sixth electron that constituted 
the formation of negative charge state was assumed to have been accepted 
from a nearby donor (presumed N^) [SI US]. The six-electron/negative- 
charge model has been generally accepted since its proposal [50], with one 
notable exception being Lenef and Rand [19] , who tentatively proposed a two- 
electron/neutral model, but later rejected it based upon Mita's assignment 
of NV*^ [73]. Direct evidence in support of the negative- charge model has 
been obtained through the observation of increases in positron trapping that 
are correlated with the formation of the NV~ optical band |145] . Having 
established the number of electrons associated with the NV charge states, 
the MO a'l that occurs within the diamond band gap will be ignored in the 
remainder of this review, since it is filled for all observable electronic levels 
of NV~ and NV°. As will be discussed in the next section, the observable 
electronic structures of the charge states are completely described by the 
various occupations of the MOs (ai, Cx-, Cy) within the diamond band gap by 
the four/ three (NV~/ NV*^) electrons that are in excess of the filled valence 
band. 
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Figure 12: The simple molecular model used by Loubser and van Wyk to 
interpret their EPR measurements of N V~ . Loubser and van Wyk [HI H5] 
explained the existence of an spin triplet state by the depicted occupation 
(arrows indicating electron spin- up/down) of the defect levels ( 
by six electrons. They also provided simple expressions (to right of defect 
levels) for the defect orbitals in terms of the sp^ orbitals (ci, C2, C3, n) of 
the nearest neighbour nitrogen and carbon atoms to the vacancy (depicted 
below the defect levels) by applying the linear combinations of atomic orbitals 
technique. 



Loubser and van Wyk did, however, make the incorrect conclusion that 
the triplet state was an excited metastable state that existed between 
optically active ^A and singlets. This conclusion was understandably 
based upon their observation that the EPR signal existed only under optical 
illumination and increased when the ZPL and phonon peaks were excited. 
Indeed, their incorrect conclusion persisted as a source of contention for over 
a decade and, as will be discussed in the next section, was ultimately re- 
solved by the observation of the EPR signal in the absence of illumination 
|146[ 1147] . The triplet state is now accepted as the ground state of the NV~ 
centre [HI [781 11461 1147i I148j and the dependence of the EPR signal strength 
on optical excitation is now understood as being due to the process of optical 
spin-polarisation [IS]. Notably, Loubser and van Wyk did recognise that the 
triplet spin sub-levels were being selectively populated under optical excita- 
tion as the high field lines in their EPR spectra were in emission. 

The optical band of NV° was initially only observed in cathodolumines- 
cence (note that NV~ is not generally cathodoluminescent) due to its weak 
photoabsorption signal [TB]. The NV° optical band was eventually observed 
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in photoabsorption and photoluminescence by Davies [7B], who noted the 
presence of the Jahn- Teller effect in the optically active states. Although 
Davies identified the optical band as an E-A electronic transition occurring 
at a trigonal centre, it was not until Mita [75] observed correlated changes in 
the concentrations of NV~ and NV° under heavy neutron irradiation, that it 
was asserted that the defect centre responsible for the NV" optical band was 
indeed the neutral charge state of the NV centre [7S]. Mita supported his 
assertion using the argument that damage induced by neutron irradiation 
lowers the diamond Fermi level, such that NV~ reduces in concentration 
relative to NV°. Mita's conclusion was generally correct, although it was 
subsequently concluded that Fermi level arguments do not provide correct 
predictions of the relative concentrations of charge states of defects in dia- 
mond and that the concentrations are instead governed by the microscopic 
distributions of donors (such as N^) in the proximity of each individual de- 
fect |149] . Consequently, in type lb diamond (substantial concentration), 
there is a large mean number of proximal donors and their microscopic 
distributions about the NV centres are reasonably uniform. Thus, NV~ is 
the most prevalent charge state in type lb diamond [S3]. Whereas in type Ila 
diamond (undetectable concentration), there is a small mean number of 
proximal donors and their microscopic distributions about the NV centres 
vary significantly. Therefore, there is no clearly dominant charge state and 
there is significant variation in the relative concentrations of the NV charge 
states in type Ila diamond [149] . The dependence of the charge state con- 
centrations on the microscopic distributions of donors also implies that for 
centres formed using N"*" ion implantation in type Ila diamond, the charge 
state concentrations depend sensitively on ion fluence, straggling and the 
other implantation and annealing parameters that determine the nitrogen 
distribution |117[ I126[ 1127] . The microscopic distributions of acceptors (such 
as the lattice vacancy and surface traps of nanodiamonds) have also been 
observed to affect the charge state of the NV centre [991 IM [M [TT2l[T5n] . 

Mita's identification of NV° was eventually confirmed by observations of 
photoconversion between the two charge states via the dependence of the rel- 
ative emission intensities of the charge state optical bands on the excitation 
intensity [H31 1103] and energy [521 1102] . As previously mentioned, it is now 
believed that photoconversion leads to an equilibrium of the charge state con- 
centrations under optical excitation and that the equilibrium concentrations 
are dependent on the excitation energy and intensity [120] . In spite of the 
observations of photoconversion and the dependence of the centre's charge 
state on the various formation factors discussed, there does not presently ex- 
ist a model that accounts for both the microscopic origins of the NV centre 
charge state and of the photoconversion process. This outstanding issue will 



20 



be discussed further at a later stage in the review. 

Before preceding to a review of the known aspects of the electronic struc- 
tures of NV~ and NV°, it is worth mentioning that a third optical ZPL in 
the range 4.325-4.328 eV has been attributed to the NV centre |151l 1152] . 
This third optical ZPL has been positively correlated with the growth of the 
NV~ optical ZPL following radiation damage and annealing in type lb dia- 
mond [1511 1152] . This third optical ZPL has not been studied in any detail, 
however it has been speculated that it corresponds to the transition of an 
electron between the valence band and the partially filled E defect levels {cx, 
Cy) |153j . If this speculation is correct, then the third ZPL will be another 
manifestation of the microscopic origins of the NV charge state and may be 
involved in the photoconversion process. 

3 Electronic structure 

The electronic structures of NV~ and NV° have been under constant debate 
since their first reports in the 1970s [TTl ITS] and several aspects are still un- 
der contention. The currently accepted electronic structures are depicted in 
figure [131 As discussed in the previous section, the orbital symmetries of the 
electronic states involved in the optical transitions of NV~ and NV" were 
established very early by photoluminescence polarisation and uniaxial stress 
studies [TTl [751 [139] . However, since the NV~ EPR signal was initially only 
detected under optical illumination [HI [15] , it took some time for the spin 
multiplicities of the NV~ ground and optically excited states to be correctly 
identified. The first evidence of the spin triplet nature of the NV~ ground 
state was obtained by Reddy et al ^8] in their optical hole burning and mag- 
netic circular dichroism study. Further evidence of a triplet ground state was 
obtained using ODMR [TS], spin-locking [7S] and spin cross- relaxation |148j 
studies. Confirmation was ultimately provided by observations of the EPR 
signal in the dark |146[ 1147] . Similar spin identification problems to NV~ 
existed in establishing the spin multiplicities of the NV° ground and optically 
excited states. The long absence of any EPR signal that was correlated to 
the NV° optical band, facilitated only speculations that the ground state was 
a spin doublet whose EPR signal was undetectable due the Jahn- Teller effect 
[86], which had been observed very early in studies of the optical transition 
|76j . The collection of sound evidence for the doublet nature of the NV° 
ground and optically excited states required the correct identification of the 
charge state [7S| and the subsequent observation of a metastable spin quartet 
EPR signal that was consistent with the neutral charge state identification 
|86j . Notably, the exact symmetry (possibly "^Ai or ^^2) of the optically ex- 
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cited doublet has not yet been determined, which is particularly important 
to the establishment of the NV'' electronic structure, since the molecular 
model predicts three doublet states {^Ai, ^^2, '^E) associated with the aie^ 
configuration. 
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Figure 13: Schematics of the known aspects of the (a) NV~ and (b) NV° 
electronic structures [16], [86]. The energies of the observed ZPL transitions 
between the electronic levels indicated. The electron occupations of 

the NV MOs associated with each electronic structure level are denoted in 
brackets. 



3.1 The fine structure of the NV level 

The ~2.88 GHz fine structure splitting of the NV" ground state was 
measured in the first EPR study (SI [15]. However, it was not possible to 
energetically order the = 0, ±1 spin sub-levels until Raman-Heterodyne 
electron nuclear double resonance (ENDOR) and double nuclear magnetic 
resonance (double NMR) measurements were performed pT| [M]. There was 
initially some debate over whether the fine structure splitting was due to elec- 
tronic spin-spin interaction or spin-orbit interaction or both. Simple symme- 
try considerations [19] imply that spin-spin interaction induces a first order 
energy splitting of ^^2, whereas spin-orbit interaction only induces a second 
order energy splitting. The observations of GHz scale spin-orbit interactions 
in the fine structure [8^ [85] . combined with the absence of any indica- 
tions of an excited state within a similar energy scale to the ground state. 
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have led to the consensus that second order spin-orbit terms are hkely to 
be too small to be principally responsible for the ground state fine structure 
splitting. Thus, the splitting has been attributed to first order electronic 
spin-spin interaction [HI HHl |5T] . 
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Figure 14: Fine and hj^erfine structures of the ^^NV~ ground state [136] . 
The fine structure levels are denoted by their spin-orbit symmetry and = 
0, ±1 spin projections and the hyperfine structure levels are denoted by their 
hyperfine symmetry and mj = 0, ±1 spin projections. 



The first EPR, ENDOR and NMR studies of NV" also provided mea- 
surements of the ground state hyperfine structure due to the 1 = 1 nuclear 
spin of the centre's native ^^N nucleus. The fine and hyperfine structures of 
the ^'^NV" ground state are depicted in figure [Ml The fabrication of NV~ 
centres using ^^N+ ion implantation has facilitated the observation of the dif- 
ferent ground state hyperfine structure arising from the I = \ nuclear spin 
of ^^N (refer to figure [T^ . The ground state fine and hyperfine structures for 
both nitrogen isotopes can be described by the canonical spin-Hamiltonian 
of trigonal defects |143] 



H, 



D 



gs 



-P, 



gs 



Sl-S{S+l)/'i\+^lsJ, + A^^ 
i!-I{I + l)/3 



So-Jt. 



Syly 



where Dgs ~ 2.88 GHz is the fine structure splitting, Pgs is the nuclear electric 
quadrupole parameter and and Aj^^ are the axial and non-axial magnetic 
hyperfine parameters. Note that since ^^N is a spin-half nucleus, it does not 
have an electric quadrupole moment. Whilst the magnitudes and signs of Dgg 
and Pgs were accepted from their first measurement H5| 1144] , the signs 
of the magnetic hyperfine parameters have been subject to debate. Table 
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[T] contains the various measured and calculated values of the NV~ ground 
state hyperfine parameters. Note that at room temperature the Dg^ reduces 
to ~2.87 GHz [1]. The temperature dependence of Dgs will be discussed later 
in this review. 
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Figure 15: Fine and hj^erfine structures of the ^^NV~ ground state [124] . 
The fine structure levels are denoted by their spin-orbit symmetry and irtg = 
0, ±1 spin projections and the hyperfine structure levels are denoted by their 
mi = ±1 spin projections. 



The debate over the signs of the magnetic hyperfine parameters has its 
origins in the unpaired spin density distribution of the NV~ ground state 
and the relative magnitudes of the Fermi contact and dipolar contributions 
to the parameters. The magnetic hyperfine parameters are expressed in 
terms of the Fermi contact (fgs) and dipolar (dgs) contributions by = 
fgs + 2dgs and A^^ = fgs — dgg |143j . The Fermi contact contribution is 
directly proportional to the unpaired spin density at the N nucleus and the 
dipolar contribution is related to the p atomic orbital character of the spin 
density in the region centred on the N nucleus [HI HS]. Using the simple 
molecular model of Loubser and van Wyk [HI US], the unpaired spin density 
is described by the occupation of the E MOs {e^, Cy) by two electrons. Since 
these MOs are not symmetric, their amplitude along the [111] symmetry 
axis of the NV centre, where the N nucleus resides, is zero. Furthermore, 
as the E MOs have no contribution from the symmetric N dangling sp^ 
orbital, the unpaired spin density produced by the E MOs is small in the 
region of the N nucleus as it is mainly located in the vicinity of the vacancy 
and its nearest neighbour carbon atoms. Therefore, the simple molecular 
model of Loubser and van Wyk suggests that the Fermi contact and dipolar 
contributions should be zero and very small, respectively, in contradiction to 
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Table 1: Experimental measurements and ab initio calculations of the hyper- 
fine parameters of the NV~ ground state tabulated by reference and nitrogen 
isotope. Uncertainty ranges are included if they were provided by the source 
references. Note that the ^^N nucleus does not have a quadrupole moment. 



Reference 


Iso. 


Ajs (MHz) 


Al (MHz) 


Pgs (MHz) 


Loubser^ [H iS] 


14N 




±2.32 ±0.01 




He'^ |144j 


14N 


+2.10 ±0.10 


+2.30 ±0.02 


-5.04 ±0.05 


Felton<= |SB] 


14N 


-2.70 ±0.07 


-2.14 ±0.07 


-5.01 ±0.06 


Steiner^ [98] 


14N 




-2.166 + 0.01 


-4.945 ± 0.01 


Smeltzer*^ |154j 


14N 




-2.162 ±0.002 


-4.945 ± 0.005 


Gah" [59] 


14N 


-1.7 


-1.7 




Rabeau'^ [l24j 


15N 


-3.1 


-3.1 


n/a 


Fuchs*^ [ig 


15N 


3.01±0.05 


3.01±0.05 


n/a 


Felton<= [130] 


15N 


3.65 ±0.03 


3.03 ±0.03 


n/a 


Gali'^ [63] 




2.7 


2.3 


n/a 


'^Experiment, "'ab initio 



observation. Given that the unpaired spin density of the valence electrons 
does not adequately account for the hyperfine parameters, the parameters 
must arise from modifications of the orbitals of the N core electrons [SH|- If 
the N core orbitals are slightly different for different electron spin-projections 
(known as core polarisation), then a significant Fermi contact contribution 
will result due to the significant amplitudes of the core orbitals at the N 
nucleus. Furthermore, the Fermi contact contribution of the core electrons 
will be negative, as opposed to the (by definition) positive contribution of 
the valence electrons |143] . Since the N core orbitals are expected to be 
principally s atomic orbital in character, the dipolar contribution of the core 
electrons is expected to be very small. Thus, the hyperfine parameters are 
dominated by the negative Fermi contact contribution of the core electrons. 
This conclusion was overlooked in the first measurements of the hyperfine 
parameters [HI EH EH 1124] and was eventually identified in an ab initio 
study ^D]. Since the conclusion was overlooked in the first measurements, 
the signs of the magnetic hyperfine parameters were incorrectly assigned as 
positive in order to be consistent with their assumed valence spin density 
origin. Later hyperfine measurements [86l [79l [98l I130[ 1154] have had the 
benefit of the identification of core polarisation and have correctly assigned 
the signs of the hyperfine parameters. 

The influence of static electric E, magnetic B and strain 6 fields on the 
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Table 2: Experimental measurements of the electronic g- factor and electric 
dipole components of the NV~ ground state tabulated by reference. 



Reference 








4s 








(Hz cm/V) 


(Hz cm/V) 


Loubser [HI gS] 


2.0028±0.0003 


2.0028±0.0003 






Felton |130| 


2.0031±0.0002 


2.0029±0.0002 






van Oort [SI] 






17±3 


0.35±0.02 



NV ground state has been observed to be well described by the addition of 
the following potential to the spin-Hamiltonian ([1]) 

Vgs = ^bqIsSzBz + HB9gs {SxBx + SyBy) + I2n9nI- B 

+dl{E, + 6z) [Sl - S{S + l)/3] + dj^iE, + 6,) (^J - Si 



+ dj^{Ey + 6y) [SxSy + SySx) (2) 

where /x_b is the Bohr magneton, ^jy is the nuclear magneton, (yf|^ and g^^ 
are the components of the ground state electronic g-factor tensor, Qn is the 
isotropic nuclear g-factor of ^"^N or ^^N as required, and and dj-^ are 
the components of the ground state electric dipole moment. The cartesian 
components of the fields are defined with respect to the coordinate system 
depicted in figure HSl Note that crystal strain has been expressed as an 
effective electric field in the above [16] and that all fields are assumed ap- 
proximately constant over the dimensions of the NV centre. Measured values 
for the ground state electronic g-factor and electric dipole components are 
contained in table [2J Being of A2 orbital symmetry, the ground state can- 
not possess a first order orbital magnetic moment or a non-axial permanent 
electric dipole moment. Consequently, the ground state electronic g-factor 
components can only differ from the isotropic free electron value Qe ~ 2.0023 
if higher order terms in the centre's interactions with magnetic fields are 
of adequate magnitude or if excited electronic states with orbital magnetic 
moments are significantly coupled (via electronic spin-orbit and spin-spin 
interactions) with the ground state |143j . Likewise, the ground state may 
only possess a permanent electric dipole moment if excited electronic states 
with electric dipole moments are significantly coupled (also via electronic 
spin-orbit and spin-spin interactions) with the ground state [2^ I136j . 
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Figure 16: Schematic of the NV centre depicting the nearest neighbour ni- 
trogen (brown) and carbon (gray) atoms to the vacancy (transparent) as well 
as the adopted coordinate system (xyz). The coordinate system is defined 
such that the z axis is coaxial with the trigonal symmetry axis of the centre, 
the origin resides at the vacancy site (the centre of the defect), and the x 
axis passes below one of the nearest neighbour carbon atoms. 

3.2 The low temperature fine structure of the NV 
level 

The low temperature fine structure of the NV~ excited state (refer to 
figure [T7|) was only recently observed directly [811 ESj. The direct obser- 
vations were made using photoexcitation spectroscopy of spectrally stable 
single centres that exhibited narrow optical linewidths. Direct observation of 
the fine structure proved elusive in earlier spectroscopic attempts due to the 
confounding effects of strain broadening in ensembles, spectral diffusion, spin 
mixing in the strain-split fine structure and optical spin-polarisation. A num- 
ber of pioneering works employed spectral hole burning [HI I155[ I157[ 1158] , 
narrow band ODMR [H2] and four- wave mixing techniques |156[ 1159] to probe 
the strain broadened optical ZPL of NV~ ensembles. Some of these works 
|155[ I156[ I157[ I158[ 1159] gave aspects of the fine structure that were eventu- 
ally observed directly. Notably, qualitatively correct models of the fine 
structure existed very early [IHl |19], but without the direct observation of 
the structure, it was not possible to correctly assign the measured features 
to the splittings predicted by the models. 

The first low temperature spectroscopy studies of single NV~ centres 
[21 11011 1116] were also unable to directly observe the fine structure of ^E 
because the centres investigated had broad optical ZPLs due to poor spec- 
tral stability. The single centres selected for the studies also typically had 
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Figure 17: (a) Low temperature zero field fine structure of the NV~ 
excited state |46]. The fine structure sphttings are indicated in terms of 
the parameters of the effective Hamiltonian Spin-orbit symmetries and 
rris = 0, ±1 spin-projections of the fine structure states are as labelled on the 
right hand side, (b) The dependence of the fine structure on the transverse 
strain energy djrsS± = dj:g.J'^l+^ as modeled by ([3]) and (jl]). The model 
parameters used to produce the plot were those of Batalov et al [H] (refer 
to table E]). 



low strain |116] . which due to spin mixing in the lower fine structure branch, 
ensured that the lower branch did not have adequate fluorescence to be vis- 
ible [85]. The combination of optical excitation and microwaves tuned to 
the ground state magnetic resonance, enabled the sufficient increase of the 
fluorescence of the lower branch fine structure levels for them to be observed 
in the later single centre studies [SU ES] . This increase in the fluorescence of 
the lower branch occurred due to the microwaves maintaining steady state 
population in all spin-projections. The poor spectral stability of the first sin- 
gle centres to be studied is believed to be due to their formation in type lb 
diamond, which contains significant concentrations of electron donors (such 
as Ng). The donors undergo photoionisation during the optical excitation 
of NV~ and produce fluctuating electric fields in the vicinity of the NV~ 
centre that broaden the optical ZPL [lOH 1116] . The spectrally stable centres 
used to directly observe the fine structure were consequently formed in type 
Ila diamond, which has very small concentrations of donors. Thus, it can 
be seen that the achievement of direct observation of the NV~ optical fine 
structure required several advances in the understanding of NV" as well as 
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Table 3: Experimental measurements of the parameters of the low tempera- 
ture NV^ fine structure tabulated by reference. Uncertainty ranges were 
not provided by the source references. 



Reference 


Di (GHz) 


Di (GHz) 


AL (GHz) 


Xi (GHz) 


Tamarat [85] 


1 


0.5 


4.4 


0.2/V2 


Batalov [H] 


1.42 


1.55/2 


5.3 


0.2/V2 



spectroscopy and fabrication techniques. 

Given the assignment of the energy ordering of the ground state spin 
sub-levels as well as the direct observation of the strain and electric field 
dependence of the low temperature fine structure [MII85] . it was straight 
forward to assign the spin-projections of the ^E fine structure levels using 
optical spin selection rules (refer to figured?]). The double degeneracy of 
the lower two levels imply that they are of E spin-orbit symmetry, which 
is in keeping with the predictions of the molecular model [l6l HTJ SHI HHl 
[52| [53] . Whilst the molecular model predicted the Ai and A2 spin-orbit 
symmetries of the remaining non-degenerate levels, their energy ordering 
had to be interpreted from experiment. The ordering of the upper levels has 
very recently been assigned as depicted in figure [T7] using lifetime arguments 
|10j . These lifetime arguments involve non-radiative processes and will be 
discussed further at a later stage in this review. 

The low temperature ^E fine structure can be described by the following 
effective Hamiltonian 



H. 



LT 



S^^ - S{S + l)/3\ -Xl&y^S, 
a, ® {Si - Si) - a., ® {SyS^ + S^Sy] 
d-z ® {S^Sz + SzS^) -cTx® {SySz + SzSy] 



(3) 



where a^, d'y and az are the standard Pauli matrices that represent the fine 
structure orbital operators in the basis of orbital states {fE^), \^Ey)} asso- 
ciated with the ^E level, and S^, Sy and are the S = 1 spin operators. 
The values of the fine structure parameters DI^, D^,, A|!^ and X^ as obtained 
from the direct observations [HH [85] of the low temperature fine structure 
are tabulated in table [31 

Applications of the molecular model ^[ [S31 [HI] clearly indicate that the 
physical origins of D^^ and D^. are spin-spin interactions and the origin of 
Aj!^ is axial spin-orbit interaction. The origin of A^^ has, however, experi- 
enced a degree of contention. Manson et al jUJ [52] asserted that A^^ arises 
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from non-axial spin-orbit interaction, whereas Lenef and Rand asserted that 
it could not arise from spin-orbit interaction using arguments based upon 
the Hermitian properties of the spin-orbit operator. It was later shown sep- 
arately by Doherty et al [l6] and Maze et al [53] that X-^^ has contributions 
from spin-spin interaction. Furthermore, Doherty et al extended Lenef and 
Rand's argument to definitively prove that spin-orbit interaction does not 
contribute to X-^^. However, the contention has not yet been completely re- 
solved because two-electron spin-orbit interactions (spin-own-orbit and spin- 
other-orbit) have not yet been considered in the molecular model and are 
likely to also contribute to Af^. 

The influence of static electric E, magnetic B and strain 6 fields on the 
low temperature fine structure has been observed to be well described by 
the addition of the following potential to the effective Hamiltonian ([3D 

Kf = dliE^ + 6.) + diiE, + 6^,)a^-diiEy + 6y)a, 

+Mt^y + 9lsSz)B, + fiBQi {SxB, + SyBy) (4) 

where l\g is the ^E orbital magnetic moment, g\g and g^r^ are the components 
of the ^E electronic g-factor tensor, and and dj:^ are the components 
of the ^E electronic electric dipole moment. Out of all of the interaction 
parameters, only /||^ has been measured directly at low temperature. The 
magnetic circular dichroism measurements of Reddy et al [IH] and Rogers et 
al [87] both obtained ll, = 0.1 ± 0.01 at ~1.5 K. The low temperature ^E 
Stark effect measurements of Tamarat et al [851 1114j indicate that the ^E 
electric dipole moment components are of the order ~6 kHz m/V. However, 
as the direction of the electric field was not known in their experiments, 
the individual components of the ^E electric dipole moment could not be 
distinguished. The components of the ^E electronic g-factor tensor have 
only been measured at room temperature and were found to be g^^ ~ gl^ ~ 
g^^ = 2.01 ± 0.01 [791 [80]. Due to the presence of the Jahn-Teller effect in 
^E, it is likely that the low temperature g-factor tensor and electric dipole 
components differ from those measured at room temperature [SB] • A detailed 
analysis of the influence of the Jahn-Teller effect on the ^E orbital magnetic 
moment, electric dipole moment and g-factor tensor is yet to be presented in 
the literature. 

3.3 The room temperature fine structure of the NV 
level 

The first observation of the room temperature ^E fine structure was made 
by Epstein et al |160j in their ODMR study of single centres and ensembles. 
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Using precise magnetic field alignment, Epstein et al detected a LAC at an 
axial magnetic field of ~500 G, which they assigned to ^E. This assignment 
was later confirmed by direct ODMR detection of the room temperature 
fine structure using single centres [791 180] . These direct observations were 
made at similar times to the first direct observations of the low temperature 
fine structure [SH However, the observations P, [70] incorrectly iden- 
tified the room temperature fine structure as the upper branch of the low 
temperature fine structure. The misidentification was justifiably made for 
two reasons. The first, was that there were only three fine structure lev- 
els observed at room temperature, which is consistent with one of the low 
temperature branches. The second reason was that the upper branch was 
known to undergo bright spin-conserving transitions in contrast to the lower 
branch. The lower branch was thus concluded to be unobservable at room 
temperature due to the presence of spin-flip transitions and increased non- 
radiative decay at room temperature. The room temperature fine structure 
was correctly identified as being the result of phonon mediated averaging of 
both low temperature fine structure branches by Rogers et al |87j . 

Rogers et al observed the temperature dependence of their CW ODMR 
detection of the excited state LAC in a NV~ ensemble. For temperatures 
down to 25 K, there existed a single LAC at ~ 510 G that was consistent 
with the room temperature fine structure observed in the single centre stud- 
ies j7^ [50] . Below 25 K, Rogers et al observed an additional LAC feature 
at ~ 250 G and that this feature varied with applied stress. By modeling 
the strain dependence of the low temperature fine structure, Rogers et al 
showed that the behaviour of the LAC features at ~ 250 G and ~ 510 G 
in response to applied stress was consistent with the low temperature fine 
structure. Hence, Rogers et al had demonstrated the transition between the 
low and room temperature fine structures. In order to explain their ob- 
servation, Rogers et al proposed that phonon transitions between the low 
temperature fine structure levels resulted in an average fine structure that 
varied with temperature, such that at low temperature the averaging effect 
was minimal and at room temperature the averaging resulted in the observed 
room temperature fine structure. Specifically, the averaging process is be- 
lieved to occur via the redistribution of population optically excited to a given 
^E low temperature fine structure level to all the other fine structure levels 
with the same spin-projection through spin-conserving phonon transitions. 
The degree of population redistribution that occurs within the lifetime of 
^E depends on the rates of the phonon transitions, which in turn depend on 
temperature. At sufficiently high temperature, the phonon transition rates 
will greatly exceed the ^E decay rate and the population will be equally dis- 
tributed between the low temperature fine structure levels of the occupied 
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spin-projection. Therefore, at sufficiently high temperatures the energy of 
the centre's state is the average of the low temperature fine structure lev- 
els of the occupied spin-projection. The equivalence of the average of the 
low temperature fine structure and the room temperature fine structure was 
con&med by the single centre experiments of Batalov et al [H] . A detailed 
model of the averaging process is yet to be presented and a thorough anal- 
ysis of the transition between the low and room temperature fine structures 
remains an outstanding issue to be resolved. 

The observed room temperature fine and hyperfine structures can be 
described by a spin-Hamiltonian that is equivalent to that used to describe 
the ground state 



^i!^ = D 



es 



~ ^{^ + + ^es^zlz + ^es S^Ix + Syly 



+Pe4/:-/(/+l)/3j (5) 

where the fine and hyperfine parameters are tabulated in table HI Note that 
the hyperfine structure has only been observed at room temperature and no 
model has yet been developed for the low temperature hyperfine structure. 
Comparing the above with H^J' ([3]), it is evident that the averaging processes 
can be described by a partial trace over the orbital states H^^'' = jtY^{H^^}, 
which is conceptually consistent. 

The influence of static electric E, magnetic B and strain 6 fields on the 
room temperature ^E fine structure has been observed to be well described 
by the addition of the following potential to the above spin-Hamiltonian (jSj 

Kf = dl{E, + 6^)+g^^S-B + aS'y-Sl) (6) 

where ^ is a parameter that has been related to strain [701 EDI Ell ; 
previously discussed, g^^ = 2.01 ±01 [791 180] is the measured isotropic room 
temperature electronic g-factor of ^E. Comparing the above with the low 
temperature potential (jlj) and extending the notion of the room temperature 
averaging being consistent with a partial trace over the orbital states, it is 
not clear where the strain related term originates from, since all of the 
non-axial strain related terms in (j4]) are traceless. Fuchs et al [79] measured 
^ ~ 70 ± 30 MHz for one of the centres they investigated and related the 
parameter to strain because the strain was very large in that particular centre 
(as indicated by a ground state strain splitting of ~6 MHz). Neumann et 
al [SU] only observed ^ > for < 20% of the centres they surveyed in bulk 
diamond. However, Neumann et al did observe a splitting consistent with 
the term in a single centre within a nanodiamond. As internal strains are 
typically large in nanodiamonds, this observation is consistent with Fuchs 
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Table 4: Experimental measurements and ab initio calculations of the room 
temperature NV~ fine and hyperfine structure parameters tabulated by 
reference and nitrogen isotope. Measurement uncertainty ranges are indi- 
cated if provided by source reference. Note that the ^^N nucleus does not 
have a quadrupole moment. 



Reference 


Iso. 


Dl (MHz) 


Ai (MHz) 


Al (MHz) 


p 

^ es 


Steiner^ [SB] 


14N 


1420 


(±)40 


(±)40 




Fuchs" [79] 


15N 


1425±3 


(±)61±6 


(±)61±6 


n/a 


Neumann*^ [80] 


15N 


1423±10 






n/a 


Gali'^ [63] 


15N 




-39.2 


-57.8 


n/a 



"^Experiment, ""ah initio 



et al's conclusion that ^ is related to strain. Consequently, it appears that 
the averaging process is not complete for centres with very large strains, 
leading to an additional splitting being observed in their room temperature 
fine structures that is otherwise absent in centres with smaller strains. Hence, 
the objective of any attempt at a detailed model of the averaging process 
should be the description of the conditions required to observe the strain 
related splitting of the room temperature '^E fine structure. 

Comparing the room temperature ^E magnetic hyperfine parameters to 
those of the ^A2 ground state (refer to table [1]), it is evident that the magnetic 
hyperfine interaction is much larger in ^E. This difference can be explained 
by once again returning to the simple molecular model of Loubser and van 
Wyk [^ US]. The unpaired spin density of ^E arises from the unpaired 
electrons occupying both the Ai (ai) and E {e^, Cy) MOs in the first excited 
configuration aic^ . Unlike the E MOs, the Ai MO has a contribution from 
the dangling N sp^ atomic orbital and thus contributes to both the Fermi 
contact and dipolar components of the magnetic hyperfine interaction [79] . 
Furthermore, since these contributions arise from unpaired valence electrons 
and not core polarisation, the contributions will be opposite in sign to the 
core polarisation contributions in the ground state [63]. Thus, the magnetic 
hyperfine parameters are expected to change sign, as supported by the ah 
initio calculation of Gali |63j . 

3.4 The fine structure of the NV^^ level 

At this stage, it will be instructive to compare the fine and hyperfine param- 
eters of NV° with those of NV" ', since they are expected to have common 
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Table 5: Experimental measurements and ab initio calculations of the NV° 
fine and hyperfine structure parameters tabulated by reference and nitro- 
gen isotope. Note that the ^^N nucleus does not have a quadrupole moment. 



Reference 


Iso. 


Diy, (MHz) 


Aj,yo (MHz) 


A^yo (MHz) 


Pnvo 


Felton*^ [SS] 


15N 


1685±5 


-23.8±0.3 


-35.7±0.3 


n/a 


Gali'^ [7] 


15N 




-23.4 


-39.0 


n/a 



"^Experiment, "afe initio 



origins. By noting the change to S* = | operators, the metastable excited 
state of NV° can be described by an equivalent spin-Hamiltonian to that of 
the ^A2 ground state (and the excited state at room temperature) of NV~ 
([T]). The NV" parameters as measured by Felton et al [SB] and calculated 
by Gali [7j are contained in table |3 It is clear that the NV° '^A2 hyperfine 
parameters are similar to the room temperature NV~ parameters. This 
can be explained by the fact that the Ai MO is unpaired for both "^^2 (ctie^) 
and ^E (aie^), which implies that similar unpaired spin density will occur at 
the N nucleus for both ^A2 and ^E due to the Ai MO being non-zero at the 
N nucleus in contrast to the E MOs. There also appears to be a strong rela- 
tionship between the electronic g-factor components of NV'' ^^2 (measured 
to be g^yo = 2.0035(2) and ^J^^o = 2.0029(2) gS]) and NV" M2 (refer to 
table [2]). This is likely due to the common orbital symmetry of the two states 
leading to similar origins of the shifts of the components away from the free 
electron value, but no analysis of the shifts has been performed to confirm. 

3.5 The NV and ^Ai levels 

Returning to NV~, the existence and energy ordering of the intermediate 
singlet states {^Ai and ^E) that are predicted by Loubser's and van Wyk's 
simple molecular model [HI US] has been a constant source of contention 
throughout the study of NV~. The early contention was fueled by the in- 
terpretation of very different storage phenomena in the centre's optical dy- 
namics as evidence of a metastable state. In their four wave mixing study, 
Redman et al [147] attributed a very long storage of ~ 0.3 s to a metastable 
state in the optical cycle that they speculated to be the ^Ai singlet (ignoring 
the ^E singlet). However, their observation is inconsistent with more recent 
measurements of the optical cycle and is likely to be related to spin-lattice 
relaxation and not a metastable state |161] . Another interpretation of a 
metastable state in the centre's optical dynamics was based upon observa- 



34 



tions of fluorescence intensity fluctuations in the first single centre study pQ. 
The presence of a metastable state was also suggested by the observations 
of 'bunching' in the fluorescence autocorrelation functions of single centres 
(refer to figure [7]) in the second single centre study [2]. 

The models of the centre's optical dynamics that were employed in the 
early single centre studies only contained three levels (corresponding to the 
and states and a metastable state assumed to be ^Ai) [21 El IH 11161 
I162j . In order to explain observations of biexponential radiative decay and 
temperature dependent fluorescence saturation, the three level models re- 
quired the metastable state to repopulate with a rate that has an Ar- 
rhenius temperature dependence and a 37 meV activation energy [2]. The 
energy of the metastable state (assumed Mi) was consequently concluded to 
be within 37 meV of ^E [2]. The subsequent use of five [ITJ US] and seven 
|163] level models nullified these conclusions regarding the metastable state 
because these more sophisticated models could explain the observations (and 
a number of other aspects of the centre's fluorescence) by correctly account- 
ing for spin-dependent fluorescence and temperature dependent decay from 
the metastable state to the ground state. 

In order to describe the intensity dependence of the observed fluorescence 
autocorrelation functions, the three level models also required the metastable 
state to optically repopulate the ^E with a rate that is linear in optical in- 
tensity |162] . This is now understood as another demonstration of the in- 
adequacy of the three level model to describe the centre's optical dynamics 
and not a property of an intermediate singlet state of NV~. The intensity 
dependence of the fluorescence autocorrelation function is intead related to 
the photoionisation process to NV° [IZl [52]. Indeed, this is an example of 
how the assignment of the metastable state to one of the singlet states was 
confounded by the presence of two storage processes in the centre's dynamics 
[m [52] : photoionisation to NV° that results in long lived storage, and non- 
radiative decay via the intermediate singlet states that results in relatively 
short lived storage. Consequently, there have been a great variety of lifetimes 
associated with the metastable state under different optical excitation condi- 
tions that are seemingly at odds with the identification of one of the singlet 
states as the metastable state. The models of the centre's optical dynamics 
and the deconvolution of the two storage processes will be discussed further 
at a later stage of the review. 

The existence of the intermediate singlet states (^Ai and ^E) was un- 
ambiguously confirmed for the first time by Rogers et al [77] through their 
observation of the 1.190 eV infrared band under optical excitation (see fig- 
ure [18]). Rogers et al used uniaxial stress to demonstrate that the infrared 
transition occurred between electronic states of E and Ai orbital symmetry. 
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Figure 18: The NV~ infrared luminescence band at two different tempera- 
tures, including its ZPL at 1.190 eV (1043 nm) and vibronic band extending 
towards lower energy, as measured using a high density NV~ ensemble [77]. 
Note that the liquid He trace has been divided by a factor of five in order to 
compare directly to the room temperature trace. The tail of the NV~ optical 
luminescence band is evident underneath the infrared band. 



Evidence that the transition occurred between spin singlets was obtained 
through the absence of any Zeeman splitting or broadening of the infrared 
ZPL. However, in the presence of a static magnetic field, the infrared fluo- 
rescence increased in correlation with a decrease in the optical fluorescence. 
These changes in fluorescence are consistent with the increased non-radiative 
population of the upper intermediate singlet state due to the mixing of the 
triplet spin-projections, thereby providing evidence that the non-radiative 
decay suspected in the centre's optical dynamics was indeed occurring via 
the intermediate singlet states. Rogers et al established that the upper sin- 
glet state had a very short lifetime (estimated ^ 30 ns) through interpreting 
the transient response of the infrared fluorescence under optical excitation. 
Rogers et al thus concluded that the lower singlet state was the longer lived 
storage state in the non-radiative decay path. Using a model of the centre's 
optical dynamics, Rogers et al further concluded that the infrared fluores- 
cence was less than expected and speculated that the infrared emission had 
significant competition from a non-radiative decay path between the singlets. 
The measurements performed by Rogers et al were later extended by Acosta 
et al |105t 1164] . who successfully measured the lifetime of the upper singlet 
to be 0.9(5) ns and the lifetime of the lower singlet, which was much longer 
and exhibited a significant temperature dependence, ranging from 462 ns at 
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4.4 K to 142 ns at 450 K. 

In spite of the comprehensive studies of the infrared band performed by 
Rogers et al and Acosta et al, the energy ordering of the singlet states could 
not be ascertained because the band was only observed in emission. Rogers et 
al speculated that ^Ai was lower in energy than based upon the assump- 
tion that spin-polarisation occurred via the emission of only symmetric (Ai) 
phonons. This spin-polarisation mechanism requires the ^Ai to be lowest so 
that the = spin sub-level of the ^A2 ground state is preferentially popu- 
lated due to its Al spin-orbit symmetry [H] . Indeed, this is still the generally 
assumed spin-polarisation mechanism, even in light of the recent conclusive 
theoretical and experimental evidence that the ordering of the singlets is the 
opposite to that speculated by Rogers et al. The theoretical evidence for the 
now accepted energy ordering of the singlet lower than the ^Ai singlet, 
was provided by a combination of improved ab initio techniques [65l [68] and 
a detailed application of the molecular model and electrostatic identities [16]. 
The experimental evidence was provided in a follow up study performed by 
Manson et al [165] , who observed shifts of the infrared ZPL emission under 
uniaxial stress that were indicative of a LAC between vibronic levels of ^E, 
and were only consistent with ^E being the lower energy singlet. Note that 
a recent ab initio calculation by Ma et al [68] added an additional element 
to the contention surrounding the intermediate singlet states by predicting 
the presence of a third singlet {^E') between the triplet states and that the 
additional E symmetric singlet was actually the upper singlet in the infrared 
transition (i.e. predicted a ^E' — infrared transition). The prediction of 
this calculation is, however, in direct conflict with the conclusive theoretical 
and experimental evidence just discussed. 

The observation of a vibronic LAC by Manson et al is also direct evi- 
dence of the Jahn- Teller effect in ^E [165j . The Jahn- Teller effect in^E and 
the vibronic bandshape of the infrared transition is yet to be analysed in 
any detail and is not well understood. Aspects of the response of the in- 
frared transition to magnetic and strain fields are also not well understood. 
The aspects are the origins of the ^E strain splitting and the absence of a 
magnetic splitting arising from the expected orbital magnetic moment of ^E. 
The simplest version of the molecular model predicts that there should be 
no strain splitting of ^E. Doherty et al [16], however, applied the molecular 
model to a higher degree of complexity and showed that the observed strain 
splitting could arise from the exchange coupling of ^E and the higher en- 
ergy singlet ^E'. The relationships between the strain splittings of ^E and 
^E (at low temperature) and their connections with the optical and infrared 
vibronic bands have not been fully explored |165j . Continuing with the con- 
vention employed by Davies and Hamer [71] to describe the stress response 
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Table 6: The stress response parameters of the optical and infrared ZPLs of 
NV~ [165] . The Al and A2 parameters correspond to shifts of the ZPLs due 
to two different symmetric (Ai) distortions of the lattice under stress. The B 
and C parameters correspond to splittings of the ZPLs due to two different 
non-symmetric (E) distortions of the lattice under stress. Each parameter is 
in units of meV/ GPa. 



ZPL 


Al 


A2 


B 


C 


Optical 1.945 eV 


1.53 


-3.91 


0.987 


1.77 


Infrared 1.190 eV 


0.484 


-0.378 


1.22 


0.693 



of the optical ZPL, the stress parameters Al, A2, B and C of the optical and 
infrared ZPLs are contained in table O The stress parameters denote the 
shift/splitting of the respective ZPLs due to different components of applied 
stress. The molecular model also predicts that should possess an orbital 
magnetic moment equivalent to that of ^E that splits the infrared ZPL in the 
presence of a magnetic field [IS]. However, no comparable splitting has been 
observed. Given the evidence of the Jahn- Teller effect in ^E via the vibronic 
LAC |165j . the most likely explanation for the absent orbital magnetic mo- 
ment is the Jahn- Teller quenching of orbital interactions [88]. Clearly, there 
is a great deal of analysis to be performed on the infrared vibronic band and 
its interactions with applied fields. The connections between the optical and 
infrared vibronic bands and the Jahn- Teller effect will be addressed more 
thoroughly in the next section. 

The observation of the NV~ infrared band established the energy sep- 
aration of the intermediate singlet states, however, it did not provide the 
energies of the singlet states relative to the triplet states. As it currently 
stands, the singlet-triplet energy separations are experimentally unknown. 
Aspects of the centre's optical dynamics provide qualitative indications of 
the singlet energies. The rapid non-radiative depopulation of the = ±1 
sub-levels of ^E, which is of the same order of magnitude as the optical spon- 
taneous decay rate [17], suggests that ^Ai is energetically close to ^E in order 
for such rapid non-radiative transitions to occur. The comparable energy of 
^Ai to ^E implies that ^E is roughly ~0.7 eV above the ground state, which 
is itself consistent with the observed relatively long lifetime of ^E. For more 
quantitative estimates of the intermediate singlet energies, one must turn to 
ab initio calculations (refer to table [7]). However, as clearly evident in table 
[71 the ab initio calculations differ greatly in their estimates of the singlet en- 
ergies. Apart from the differences in calculation method, geometry and size. 
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the inconsistencies of the calculations have their origins in the limitations of 
current ab initio methods. 

3.6 Ab initio calculations of the NV~ electronic struc- 
ture 

The vast majority of calculations have been performed using density func- 
tional theory (DFT) or a variant. The key idea of DFT is to avoid the 
computational complexity of directly solving for the many-body electronic 
wavefunction and instead solve for the simpler ground state electronic charge 
density p(r). This idea is supported by the fact that many quantities can 
be expressed as functional of p(r). For instance, the ground state electronic 
energy of a system can be described by the functional 

E[P] = Ts[p] + / Vir-)pir)df+ [ ^P^drdr^ + ^p] (7) 

where Ts[p] is the (calculable) kinetic energy of the non-interacting system 
of electrons with the same density, V{r) is the electrostatic electron- nucleus 
attraction energy, the third term is the electrostatic electron-electron repul- 
sion energy and E^^^'^^[p] is the exact, albeit unknown, exchange-correlation 
functional that captures all of the remaining contributions to the total elec- 
tronic energy. The non-interacting electronic Kohn-Sham system |166] with 
the same density is useful as its wavefunction is simple, being typically given 
by one or a few Slater determinants constructed from a set of Kohn-Sham 
orbitals and as its kinetic energy Ts[p] captures a large part of the un- 

known kinetic energy of the interacting system. The central approximation 
that is made in the application of DFT is the selection of an approximate 
exchange-correlation functional E^^^'^°^[p], such as the Local Spin Density 
Approximation functional or the Becke-Perdew functional. Then the appli- 
cation of DFT proceeds by the minimisation of E[p] with respect to p by 
varying the orthonormal thereby resulting in the approximate total 

energy and charge density of the electronic ground state. 

Time-dependent DFT aims to go beyond ground state properties, ob- 
taining information on excited states through studying the response to a 
time- dependent potential V^(r, t). However the analogue to E'^pp™^ is more 
complex and, being a more recent development, has been less studied. Al- 
ternatively, Gunnarsson and Lundqvist |167] have shown that DFT can be 
extended beyond the overall ground state to the lowest state of each distinct 
spin and spatial symmetry, e.g. allowing estimates of the ^E, ^Ai and ^E 
states in addition to in the NV centre. However though E^dp] should 
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Table 7: Summary of ab initio calculations of the NV~ electronic structure. 
The Method column details the calculation method (DFT: Density Func- 
tional Theory, TD-DFT: Time Dependent DFT, CI: Configuration Interac- 
tion, CASSCF: Complete Active Space Self Consistent Field, DFT-MBPT: 
Combined DFT and Many Body Perturbation Theory) and exchange- 
correlation functional for DFT methods (LSDA: Linear Spin Density Ap- 
proximation, BP: Becke-Perdew, PBE: Perdew-Burke-Ernzerhof, GGA: Gen- 
eralised Gradient Approximation, HSE - Hartree-Fock hybrid). The Size 
column details the number of carbon atoms (clusters) or total number of 
atoms (supercells) in the geometries employed. All energies (except for those 
denoted) were calculated in the ground state nuclear equilibrium geometry. 



State 


Reference 


Method 


Geometry 


Size 


Energy (eV) 


11/ 




DFT T,SDA 




ou 


n AA 

U .tit: 








OUpcl Cell 




r^jC\ Q 

~u.y 






DFT-RP 




4-9 984 


n 499 489 






CI 


Cluster 


42 


0.629-0.644= 




Ma [68] 


DFT-MBPT 


Supercell 


256 


0.40 


'A, 


Goss [55] 


DFT-LSDA 


Cluster 


36 


1.67 




Gali [Sg 


DFT-LSDA 


Supercell 


512 


~0.0 




Delaney [HS] 


DFT-BP 


Cluster 


42,284 


2.096,2.028 






CI 


Cluster 


42 


2.060 




Ma [68] 


DFT-MBPT 


Supercell 


256 


0.99 




Goss [55] 


DFT-LSDA 


Cluster 


36 


1.77^^ 




Larrson [5B] 


DFT-BP 


Cluster 


163,284 


1.847,1.867 




Gah [Sg 


DFT-LSDA 


Supercell 


512 


1.91 




Hossain [62] 


DFT-GGA 


Supercell 


256 


1.912 




Lin [H] 


TD-DFT 


Cluster 


85,89,104 


2.17'' 






CI 


Cluster 


85,89,104 


2.66'' 






CASSCF 


Cluster 


85,89,104 


2.53'' 




Gali [Si] 


DFT-PBE 


Supercell 


512 


1.910 






DFT-HSE 


Supercell 


512 


2.213 




Delaney [65] 


DFT-BP 


Cluster 


42,284 


1.270,1.898 






CI 


Cluster 


42 


1.932-1.958'^ 




Ma [68] 


DFT-MBPT 


Supercell 


256 


2.32 




Larrson [SB] 


DFT-BP 


Cluster 


163,284 


2.277,2.296 




Gali [Sg 


DFT-LSDA 


Supercell 


512 


~2.41 




Ma [68] 


DFT-MBPT 


Supercell 


256 


2.25 


'^ZPL 


energy; ''avera^ 


le over cluster 


sizes; '^split level due to Cs symmetry 
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then depend on the symmetry no prescription was given for this dependence, 
and codes typically implement symmetry-independent E^clp] which can lead 
to nonsensical results (e.g. for atomic carbon |168] ). Fortunately, as explained 
by von Barth |168] . symmetry-independent E^c give reasonable results for 
states which reduce in the non-interacting system to single Slater determi- 
nants of Kohn-Sham orbitals. Now in the NV centre, the '^A2 {rris = ±1) and 

states can be approximately written as single Slater determinants, en- 
abling e.g. Goss's 1.77 eV estimate of the NV~ ^7^2 -^^E vertical excitation 
energy [501 [55] . However each of the non-interacting (m^ = 0), ^E and 
^Ai wavefunctions needs two determinants. This is a common issue whenever 
two electrons are to be distributed between two orthogonal spatial orbitals, 
here and Cy. von Barth also proposed to apply DFT to calculate the en- 
ergy expectation values of non-eigenstate single Slater determinants formed 
from well-defined mixtures of symmetries, such as -^(^ Ey+^ A2{ms = 0)) or 
■^(^Ai—^Ex). The energies of these orthogonal mixtures will give an average 
of the desired energies of the components, and given the ^Ao energy those 
of ^E and ^Ai follow, von Barth's procedure has been applied a number 
of times in calculations of the NV~ centre's electronic structure and, if the 
calculation performed by Gali et al [59] is ignored, von Barth's procedure has 
enabled the intermediate singlet states to be correctly energetically ordered 
by DFT calculations [551 165] . 

The Configuration Interaction (CI) method is a more rigorous means of 
calculating the energies of excited states than DFT. CI attempts to directly 
solve for the many-body wavefunction of each electronic state by expanding 
the state in terms of the many Slater determinants formed from an expansive 
basis of single electron orbitals. The method therefore avoids the inherent 
restrictions of DFT and can in theory calculate the energy of any excited 
state at the cost of a significant increase in the required computational re- 
sources. Consequently, CI calculations are currently limited to small cluster 
models where the confinement of the defect wavefunction potentially per- 
turbs the calculated energies. Indeed, the effects of confinement are believed 
responsible for ^Ai being predicted to have higher energy than ^E in the CI 
calculation performed by Delaney et al [55] because the ^E energy was found 
to be dependent on cluster size. The CI calculation performed by Delaney 
et al obtained the correct ordering of the intermediate singlet levels, thereby 
resolving the disagreement in the orderings estimated by applications of von 
Barth's DFT procedure. The correct ordering of the intermediate singlets 
states has also been obtained using various quantum chemical methods em- 
ployed by Zyubin et al [BU] and post-DFT methods employed by Ma et al 
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The complementary nature of the different theoretical methods and ex- 
periment is clear in how the results of the current calculations support the 
ordering the centre's electronic states inferred from the molecular model |16] 
and experiment [165] . Furthermore, ab initio calculations have provided a 
direct means to confirm the structure of the NV centre, its charge states, and 
aspects of the molecular model, such as the number and symmetry of MOs 
associated with the centre within the diamond bandgap. 

It is important to note in interpreting the calculated energies of table [7] 
that they correspond to the energies of the electronic states in the ground 
state nuclear equilibrium coordinates (i.e. vertical excitation energies) and 
are not directly related to the observed optical and infrared ZPLs. To repro- 
duce the ZPL energies, the Stokes and anti-Stokes shifts associated with the 
correlated relaxations of nuclear coordinates with electronic transitions must 
be considered. These shifts, alongside the other observed vibronic effects in 
the NV centre, will be discussed in the next section. 

4 Vibronic structure 

4.1 The Huang- Rhys model of the NV optical vi- 
bronic band 

The vibronic structure of the NV centre is generated by the interactions of 
its electronic and vibrational degrees of freedom. The observable aspects 
of vibronic structure are most evident in the coupled electronic and vibra- 
tional changes that occur in the process of an electronic transition, such as 
the optical transition of the NV~ centre. The simplest approximate descrip- 
tion of these changes for a particular electronic transition is provided by the 
Huang- Rhys model |143j . within which, the quasi-continuum of vibrational 
modes of the diamond lattice are replaced by a single effective vibrational 
mode whose normal coordinate corresponds to the displacement of the nuclei 
between their equilibrium coordinates associated with the electronic states 
of the transition. As depicted in figure dni the use of an adiabatic approx- 
imation (such as the Born-Oppenheimer approximation) results in the vi- 
brational motion associated with a given electronic state being described by 
the potential formed by the state's electronic energy and its dependence on 
the nuclear coordinates |169] . The harmonic approximation facilitates a fur- 
ther simplification through the approximation of the potential by a simple 
quadratic function of the nuclear coordinates |169] . which for the and 
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electronic states of NV take the form 

E¥i{Q) = E3E + aQ + ^{uj^ + b)Q^ 

= E3E-En + l{uj^ + b){Q-6Qy (8) 

where Ese is the electronic energy of the excited state at the ground state 
nuclear equilibrium coordinate [Q = 0), a and b are the linear and quadratic 
electron-vibration coupling parameters of the ^E state, respectively, u and 
\/uj^ + b are the frequencies of the mode Q in the and ^E states, respec- 
tively, 6Q = —a/{uj'^ + b) is the equilibrium displacement of the ^E state, 
Eji = a^/2{uj'^ + b) = Shy/u'^ + b is the relaxation energy of the ^E state and 
S is the Huang-Rhys factor of the optical transition. The physical origins 
of the linear and quadratic electron- vibration coupling terms in i?3^(Q) can 
be understood as being the nuclear displacement and vibrational frequency 
shift that result from the redistribution of electronic charge in the excited 
electronic state |169] . Given the harmonic potentials, the vibrational states 
associated with the and ^E electronic states will be simple harmonic 
oscillator states with energies huj{i' + |) and hy/co^ + b{u + |), respectively, 
where u is the vibrational occupation |140] . 

Applying the Frank-Condon principle |169] . at low temperature the elec- 
tronic transition from the ground state to the optically excited ^E state 
originates from the ground vibronic level of ^yl2, it approximately occurs 
within the nuclear equilibrium coordinates and is followed by a non- 
radiative relaxation to the ground vibronic level of ^E that effectively dis- 
places the nuclei to their ^E equilibrium coordinates. The electronic tran- 
sition from ^E back to proceeds in the opposite sense. Due to the 
difference in the nuclear equilibrium coordinates (and to a lesser degree the 
difference in the vibrational frequencies) associated with the two electronic 
states, the overlap of the vibrational wavefunctions of the two electronic 
states are non-zero for vibrational states with different vibrational occupan- 
cies |143j . Consequently, transitions are allowed between all of the vibronic 
levels of the two electronic states with rates that are proportional to the 
vibrational Frank-Condon factors, thereby producing the vibronic band of 
the electronic transition that is Stokes shifted in excitation and anti-Stokes 
shifted in decay |143] (see figure [19]). The transitions between vibronic levels 
with the same occupation of the effective mode contribute to the ZPL of the 
electronic transition. 

Ah initio studies have calculated the nuclear equilibrium coordinates of 
the ground and excited electronic states of the optical transitions of NV~ 
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Figure 19: Configuration coordinate diagram depicting the vibronic transi- 
tions of tlie Huang- Rliys model in tfie Frank-Condon picture. Ef^{Q) and 
Es^^lQ) are tlie electronic energies of the NV~ ground and optically 
excited electronic states as harmonic functions (curved lines) of the effec- 
tive nuclear displacement coordinate Q. The horizontal lines represent the 
vibrational levels of each electronic state. AS and AAS are the Stokes and 
anti-Stokes shifted energies of the vertical optical transitions (solid arrows), 
respectively. The diagonal transition (solid arrow) is a zero-vibration transi- 
tion and contributes to the optical ZPL. The non-radiative transitions within 
the vibronic structures of the electronic states that follow the optical transi- 
tions are depicted as dashed arrows. Note that only the vertical transitions 
originating from the ground vibronic levels and the zero-vibration transition 
between the ground vibronic levels have been depicted for reasons of clarity 
and that transitions between all of the vibronic levels are indeed allowed. 
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[SaESllSSlEHlEDlEaiSSlEniand NV° [3 ED]. Referring to figure [201 the 
positions of the nearest neighbour nuclei to the vacancy can be specified rel- 
ative to the positions of the next-to-nearest neighbours by the parameters 
contained in table El Note that the ab initio calculations consistently indi- 
cate that the next-to-nearest neighbour nuclei and beyond undergo negligible 
relaxation from their positions in a defect-free crystal lattice. As determined 
by the ab initio studies performed by Gali et al the ground state equilib- 
rium positions of NV~ correspond to a symmetry preserving outwards (away 
from the vacancy) relaxation of the nearest neighbour nitrogen and carbon 
nuclei from their positions in a defect-free crystal lattice. In the ground state 
relaxation, the nitrogen relaxes further than the carbon nuclei. The ex- 
cited state equilibrium positions of NV~ correspond to a counter relaxation 
of the carbon and nitrogen nuclei from their ground state positions, where 
the nitrogen moves further inwards than the carbon nuclei move outwards 
[59] . In a separate study, Gali et al [7] calculated that the ground state equi- 
librium positions of NV° corresponded to a lowering of the centre's symmetry 
from to Cih, which is indicative of the Jahn- Teller effect in the degen- 
erate "^E ground state. Ab initio calculations [SHI EH EB] have also obtained 
the energies of the NV~ electronic states in each nuclear configuration, which 
can be used to estimate the optical ZPL energy and associated Stokes and 
anti-Stokes shifts (refer to table [Hj). The calculated Stokes shifts may be used 
to relate the electronic energies in table [7] to the observed optical ZPL energy 
of NV". 

4.2 The linear symmetric mode model of the NV op- 
tical vibronic band 

To go beyond the overly simple Huang-Rhys model and attempt to quanti- 
tatively describe the vibronic bands of NV~, one must consider the quasi- 
continuum of crystal vibrational modes [143] . Generalising the expressions 
([8]) of the Huang-Rhys model to the many mode case, the harmonic nuclear 
potentials of the ground and ^E excited electronic states of NV" are 

m 

where Qj, a^, ha and Ui are the normal coordinates, '^E linear and quadratic 
electron- vibration interaction parameters and frequency of the i^^ vibrational 
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Table 8: Summary of ab initio calculations of the NV~ geometry tabulated 
by reference and electronic state. The unrelaxed geometry of the NV centre 
as calculated using the empirical diamond bond length 1.54 A is provided in 
the first line for comparison. All distances are provided in units of A and 
the angle LCj^N between the C^-N bond and the centre's symmetry axis is 
provided in units of degrees. 



State 


Reference 


Distance (A) 










Angle (°) 






NCn 


VN 


NC 


VC 


CCc 


CC'c 


LCnN 


Unrelaxed 


1.54 


1.54 


2.51 


1.54 


1.54 


1.54 


109.5 


'A, 


Goss [55] 


1.44 








1.45 


1.45" 






Luszczek [56] 


1.50 




2.72 




1.48 


1.51 






Larrson [58] 


1.48 




2.77 




1.51 


1.51 


105.1 




Gali [SH] 


1.46 


1.68 




1.62 


1.50 


1.50" 






Zyubin [5U] 


1.49 


1.68 




1.61 


1.52 


1.54 






Delaney [65] 


1.48" 




2.92" 




1.50" 


1.50" 


105.5" 




Zhang [70] 




1.70 




1.65 










Gali [59] 
Ma [BH] 
Zhang [TO] 




1.61 

-0.05^ 

1.63 




1.67 

+0.06^ 
1.68 









"did not differentiate CCc and CC'q] ^displacement from value; 
"small cluster. 



Table 9: Experimental measurements and ah initio calculations of the NV~ 
optical ZPL, Stokes and anti-Stokes shifted energies tabulated by reference. 
Note that the difference in the vibrational zero point energies of the 
and electronic states has been neglected in the calculated optical ZPL 
energies. This is equivalent to the neglect of quadratic electron-vibration 
coupling. 



Reference 


Stokes (eV) 


ZPL (eV) 


Anti-Stokes (eV) 


Davies" [71] 


2.180 


1.945 


1.760 


Gah* [Sg 


1.910 


1.706 


1.534 


Gali^ [61] 


2.213 


1.955 


1.738 


Ma'' [sg 


2.32 


2.09 


1.95 



"Experiment; ^ah initio 
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Figure 20: Schematic of the geometry of the NV centre depicting the nuclei 
that exists within the a:2;-plane of the coordinate system defined in figure dni 
and the carbon nuclei C'q that are bonded to the carbon nearest neighbour 
of the vacancy and exist out of the xz-plane and within the xy-plane. Given 
the symmetry of the centre, the depicted elements completely describe 
the geometry of the centre. 



mode, respectively, and bij i ^ j is the quadratic electron-vibration interac- 
tion parameter that couples the i^^ and j'*^ modes. Note that Q represents the 
collection of normal coordinates of all modes. The eigenmode displacement 
coordinates of a diamond crystal containing a NV centre will have definite 
C^i, symmetry due to the reduction of the cubic and translational symme- 
tries of the perfect lattice by the presence of the NV defect. Simple symmetry 
arguments imply that state may have linear and quadratic interactions 
with both symmetric and non-symmetric modes. The electron-vibration in- 
teraction of non-symmetric vibrations and electronic states is the root of the 
Jahn- Teller effect |143] . Due to the complexity of the many-mode Jahn- Teller 
effect, linear and quadratic interactions with non-symmetric modes are typ- 
ically ignored until it is apparent that symmetric modes cannot adequately 
describe a vibronic band [140] . A further complexity typically ignored is 
the quadratic coupling of Ai modes |140j . Hence, the simplest theory that 
attempts to account for the quasi-continuum of crystal vibrational modes. 
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considers just the linear interactions with Ai modes. The general theory of 
such linear symmetric mode models was developed by Maradudin |17U] and 
applied extensively by Davies |171] to describe the optical band of NV~ and 
other colour centres in diamond. Davies and Hamer [71] later extended the 
model to include their nitrogen-tunneling explanation of the slight asymme- 
try between the NV~ optical absorption and emission bandshapes (refer to 
section 12]). The simple linear symmetric mode model nevertheless provides a 
reasonable description of the optical band and also provides a wealth of in- 
formation on the electron- vibration interactions in the ■^^42 and electronic 
states. 

Using ensemble measurements of the absorption band, Davies obtained 
an expression for the function g{uj) = h'^uj'^S{uj)pAi{(^) (where S{u) is the 
average Huang- Rhys factor and pAi{u}) is the density of symmetric modes of 
frequency u) by employing a trial and repeat method, in which he initially 
defined a trial function, convolved it with the liquid-nitrogen temperature 
zero-phonon lineshape, compared the calculated band with the observed band 
and repeated the process until a good fit was obtained. The function that 
Davies obtained is contained in figure [21] and the individual components of 
the absorption sideband obtained using the function are depicted in figure [221 
The function also provided an excellent fit for the temperature dependence 
of the normalised ZPL strength via the simple expression derived in [171] . 

By noting the definition of S{uj) = a'^{u) /2Tluj'^ (where a^{(jj) is the average 
of of over all Ai modes with frequency w), it can be concluded that the 
function obtained by Davies g{(jo) = h^uj'^S{uj)pA-^{^) = ha'^{u)pA^{ui)/2u 
is the average linear interaction strength of the electronic state with 
Ai modes of frequency u. Consequently, g{u) provides direct insight into 
the linear interactions of ^E with Ai modes. The peak of the function at 
~67 meV indicates that Ai modes of these frequencies interact particularly 
strongly with ^E. Given that the electronic states of the NV centre are 
highly localised to the nitrogen and nearest neighbour carbon nuclei to the 
vacancy, the displacements of the Ai modes with energies close to ~ 67 
meV must also be localised to some degree to these nuclei, and the modes 
may be identified as pseudo-local modes. Note that true local modes that 
involve just the displacements of the nearest nuclei to the centre must have 
frequencies greater than the maximum frequency of the defect-free lattice. 
As there are no such local mode features present in g{u), Davies concluded 
that the NV centre has no true local modes |171] . As indicated in figure [JT] 
the peak of g{uj) at ~ 155 meV is consistent with the peak in the diamond 
density of phonon modes at a similar energy. This peak in the density of 
modes corresponds to the optical modes of the lattice. Thus, the ~ 155 meV 
peak can be attributed to the increase in the density of Ai modes and is not 
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Figure 21: The symmetric mode model of the NV~ optical (a) absorption 
band and (b) normalised ZPL intensity change with temperature |171j . The 
dashed ( ) curve in (a) indicates the absorption baseline. Note the iden- 
tified 594 nm feature is the ZPL of another defect and is unrelated to the 
NV centre, (c) The linear electron-vibration interaction strength with Ai 
modes g{uj) (solid curve) as obtained from the absorption band using the lin- 
ear symmetric mode model |171j and the diamond density of phonon modes 
up{uj) (dashed curve) [172] in normalised units. The vertical dotted arrow 
marks the peak in g{u}) at ~ 67 meV and the correlated feature marking the 
change from the Debye functional form oc of the density of modes. 
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necessarily indicative of any significant localisation of the optical modes. 




AE (meV) 

Figure 22: The normalized intensities /^(Ai?) of the n = 1,2. . .8 phonon 
sidebands (solid curves) and total absorption vibronic band I{AE) = 
Y^^In^AE) (dashed curve) of NV~ as calculated using the function gioj) 
obtained by Davies |171j using the symmetric model model (refer to figure 
I^T]) . AE denotes the energy difference from the optical ZPL of NV~ and 
the sidebands begin with the n = 1 sideband on the left hand side and in- 
crease sequentially towards the right hand side. The integrated intensity or 
the v}^ phonon sideband is normalized to e~^S^/n\, such that the integrated 
intensity of the total vibronic band is one. For further details, refer to Ref. 

A simple physical explanation for the locations of the peak interactions 
in g{u}) can be obtained by considering the phonon modes of the defect-free 
lattice. The nearest neighbour nitrogen and carbon nuclei to the vacancy 
occupy equivalent lattice sites in the diamond bravais lattice, each being con- 
nected by a primitive lattice vector. Consequently, the nearest neighbours 
only undergo relative displacements to each other for phonons with recip- 
rocal lattice vectors (fc-vectors) away from the Brillouin zone centre (the F 
point). The nearest neighbours achieve maximum relative displacement for 
phonons at the boundaries of the Brillouin zone. Indeed, phonons with k- 
vectors at certain high symmetry points on the zone boundary correspond to 
vibrations where pairs of the nearest neighbours vibrate 180° out of phase. 
Hence, the highly localised electronic states of the NV centre are expected 
to interact most strongly with phonons close to the Brillouin zone boundary. 
As indicated in figure EH the lowest energy phonons in the vicinity of the 
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zone boundary have energies close to ~ 67 meV, which is indicated by the 
change in the functional form of the density of modes from the simple Debye 
form oc bJ^ |173j . The modes that reach the zone boundary close to ~ 67 
meV are the transverse acoustic phonon modes |173] , and are thus the modes 
that mix to form the identifiable pseudo-local modes in giiS). The longitu- 
dinal acoustic and optical phonon branches approach the zone boundary at 
various energies between ~ 67 meV and ~ 160 meV |173j . Consequently, 
these modes do not produce a distinct peak in g{ixi\ but instead transform 
the sharp features present in the density of modes into the smooth curve 
of g[oj) between ~ 67 meV and ~ 155 meV. The transverse optical phonon 
modes reach the zone boundary with energies ~ 155 meV |173j . however 
given the smaller peak amplitude at ~ 155 meV compared to ~ 67 meV 
in g{ijj\ it appears that in spite of the corresponding peak in the density 
of modes, the transverse optical modes must not interact as strongly as the 
transverse acoustic modes with the electrons of the NV centre. This quali- 
tative explanation for the form of giiS) has been supported by an ah initio 
calculation [7D] that obtained the density of phonon modes for a NV~ super- 
cell model and used the displacement coordinates of the phonon modes to 
produce a projected density of modes that was weighted by the amplitudes 
of the displacements of the vacancy nearest neighbours. 

In addition to his observations of the NV~ optical absorption bandshape 
and the temperature dependence of the ZPL intensity, Davies also observed 
a temperature dependent homogeneous broadening and energy shift of the 
ZPL. These effects cannot be explained by the linear symmetric mode model 
and require the inclusion of at least quadratic interactions with symmet- 
ric modes. In order to avoid the complexity of a rigorous treatment of 
quadratic electron- vibration interaction, Davies (following Maradudin |170] ) 
assumed that the quadratic interaction parameters for modes were small 
{hij/uJiUJj <^ 1) and that they could be expressed as the product of linear in- 
teraction parameters hij oc aiUj. The first assumption is supported by recent 
ah initio studies [691 [70] that have calculated only small differences (< 2%) in 
the vibrational energies between the and electronic states. This sec- 
ond assumption can only be justified by the phenomenological evidence that 
it has provided good descriptions of the vibronic bands of other defects, since 
there is no fundamental reason for the quadratic and linear interaction pa- 
rameters to be related in such a way |170j . Using these assumptions, Davies 
derived expressions for the temperature dependent homogeneous broadening 
and energy shift of the ZPL [171j . Davies was able to obtain a good fit of the 
temperature dependence of the ZPL energy. However, due to the significant 
inhomogeneous broadening present in his ensemble sample, Davies could not 
directly apply his expression for the homogeneous width, but instead had 
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to fit a convolution of a homogeneous lineshape (with a temperature de- 
pendent width given by the above) and an inhomogeneous hneshape (with 
a width with an unknown temperature dependence) to the observed zero- 
phonon hneshape at each temperature. Although Davies obtained a good fit 
using this method, this cannot be taken as evidence that the symmetric mode 
model correctly described the broadening of the ZPL, because the temper- 
ature dependence of the inhomogeneous width could have compensated for 
any inaccuracy in the symmetric mode model during the fitting procedure. 
This point is particularly pertinent since the recent direct measurements of 
the homogeneous linewidth and its temperature dependence by Fu et al [88] 
confiict with those of Davies |171] and can only be explained by the presence 
of the Jahn- Teller effect. This and other evidence of the Jahn- Teller effect 
will be discussed further shortly. 

4.3 The asymmetry of the NV optical absorption and 
emission bandshapes 

Although Davies' symmetric mode model provided good descriptions of the 
NV~ optical absorption bandshape and the temperature dependence of the 
ZPL intensity and energy, it could not describe the asymmetry of the ab- 
sorption and luminescence bandshapes observed later by Davies and Hamer 
[71] . The slight asymmetry observed by Davies and Hamer consisted of 
a small (< 10 meV) splitting of the one-phonon peak in absorption that 
does not occur in luminescence as well as broadening of the two- and three- 
phonon features in absorption compared to luminescence (refer to figure [TT]1 . 
Given the approximations and assumptions made in the derivation of the 
symmetric mode model, the failure of the model to describe the band asym- 
metry could arise from one of the following: (1) neglect of interactions with 
non-symmetric modes; (2) neglect of quadratic interactions with symmetric 
modes in the model of the bandshapes; or, (3) the failure of the harmonic 
approximation. Davies and Hamer considered both the first and third expla- 
nations for the failure of the model. Given that recent ab initio calculations 
[69| 170] as well as Davies' models of the ZPL homogeneous broadening and 
energy shift with temperature [171] all indicate that quadratic interactions 
with symmetric modes are weak, the second explanation is the least probable 
out of all three. Davies and Hamer ultimately adopted the third explanation. 

As previously discussed, Davies and Hamer attempted to explain the 
asymmetry in the absorption and luminescence bands by the notion that 
the nitrogen nucleus tunnels between its substitutional lattice site and the 
vacancy site whilst preserving the centre's symmetry, which is analogous to 
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the inversion tunneling of nitrogen in ammonia. The notion of the tunnehng 
nitrogen nucleus is too complicated to theoretically treat if the many modes 
of the lattice are considered. Instead, Davies and Hamer had to consider 
just one symmetric mode (similar to the Huang-Rhys model detailed at the 
beginning of this section) with an effective vibrational energy of ~ 67 meV 
that describes the tunneling behaviour. Although Davies and Hamer did not 
explicitly state this at the time, the tunneling model represents the failure of 
the harmonic approximation, since the nuclear potential is no longer consid- 
ered to be a simple quadratic function of the nuclear displacement coordinate 
of the one symmetric mode, but a double-well composed of a barrier between 
two quadratic functions. As depicted in figure |23l the double-well potential 
directly yields the splitting of vibronic levels with energies that are close to 
the top of the central barrier due to the overlap of the corresponding vibra- 
tional wavefunctions of each well (i.e. inversion doubling). Consequently, it 
is straight forward to conclude that if the barrier energy in was small 
enough such that the one-phonon levels at ~ 67 meV were split due to inver- 
sion doubling, then the splitting of the one-phonon absorption peak could be 
simply explained |174] . Furthermore, the absence of any obvious splitting of 
peaks in the luminescence band would require the barrier energy in to 
be higher than vibronic levels occupied by at least six phonons with energies 
of ~ 67 meV |174] . Clearly, the tunneling model was elegant in its expla- 
nation of the vibronic splitting in the absorption band, however it required 
very specific values for the barrier energies to be quantitatively correct and 
generated implications for other aspects of the optical transition that have 
since been unobserved or demonstrated to be incorrect. 

The tunneling model proposed by Davies and Hamer was developed fur- 
ther by the much later works of Kilin et al [174] and Lin et al |175] . Kilin 
et al fitted the vibronic bands of single centres using a double-well potential 
and found that the parameters of their model varied significantly between 
the individual centres, making it improbable that any features attributable 
to the model could still be visible in the ensemble measurements of Davies 
and Hamer [52]. Lin et al refined the approach by Kilin et al and found 
that in order to match the vibronic splitting of the one-phonon peak in the 
absorption band, the ZPL also exhibited a small vibronic splitting of ~ 0.6 
meV. The studies of single centres with ZPL widths much smaller than ~ 0.6 
meV following the work of Lin et al have clearly showed that no such splitting 
exists [52]. Supporting the empirical evidence against the tunneling model, 
the early ab inito calculation of Mainwood [31] predicted tunneling barrier 
energies much greater than allowed by the model. Mainwood's calculation 
has been supported by later and more sophisticated calculations [601 ES] that 
have estimated the barrier energies associated with both and to be in 
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Figure 23: Schematic of the nuclear potentials (solid curves) of the nitrogen 
tunneling model of NV^ corresponding to the ground and optically ex- 
cited electronic states. The horizontal lines indicate vibronic levels of each 
electronic state and are labelled by their vibrational number. Vibronic levels 
marked with a prime are levels that have been split due to inversion doubling 
arising from the overlap of the vibrational wavefunctions of the left and right 
potential wells near the top of the respective barrier energies. Davies and 
Hamer [71] proposed that the barrier energy of is much smaller than ^A2, 
resulting in a vibronic splitting of the one-phonon peak (and ZPL |175j ) in 
absorption and no observable splitting in luminescence. 



excess of 4 eV, which is two orders of magnitude greater than allowed by the 
model. Hence, a body of conclusive evidence against the tunneling model has 
been collected and the model is no longer accepted as an explanation of the 
band asymmetry. In light of recent observations of the Jahn- Teller effect in 
^E [SHI I14H 1142] , the attention has now been turned to a possible Jahn- Teller 
explanation for the differences in the absorption and luminescence bands. 
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4.4 The Jahn- Teller effect and the temperature depen- 
dence of the NV optical ZPL 

Although Davies and Hamer rejected the first possibihty of a Jahn- Teller ex- 
planation, their investigation of it deserves some discussion. They conducted 
ensemble measurements of the polarisation of the luminescence band at low 
temperature (17 K) using broadband optical illumination from a high pres- 
sure mercury vapour lamp. In addition, the measurements were conducted in 
the presence of a large stress applied along the [001] direction, such that the 
two branches of the optical ZPL were split by ~ 6 meV (refer to figure [2^ . As 
expected, the lower energy ZPL branch (line b in figure is predominately 
cr-polarised (perpendicular to stress axis), whereas the higher energy ZPL 
branch (line a in figure [24|) is predominantly vr-polarised (parallel to stress 
axis). Davies and Hamer showed that the small depolarisation of the lower 
branch could be explained by the random crystal strain distribution present 
in their sample. The polarisation of the luminescence band depicted in figure 
[21] represents the polarisation of the band once corrected for the depolarisa- 
tion due to random strain inferred from the ZPL. Since the emission from 
the upper branch is negligible (< 1%) compared to the emission from the 
lower branch, the polarisation of the luminescence band is governed by the 
emission from the lower branch. Consequently, in the absence of the Jahn- 
Teller effect, where optical selection rules are degraded by vibronic couplings 
within the level arising from interactions with non-symmetric modes, the 
polarisation of the luminescence band should be the same as the lower branch 
(o"-polarisation). Davies and Hamer observed that the one-phonon peak was 
only slightly depolarised (< 10% vr-polarisation) and thus concluded that 
the Jahn- Teller effect was very weak in and not capable of explaining the 
differences in the absorption and luminescence bandshapes. 

It is an important point to note that even though the ensemble was being 
excited with broadband illumination in the experiment performed by Davies 
and Hamer, the luminescence of the upper branch is much weaker than the 
lower ZPL branch. This implies that any population excited to the upper 
branch rapidly decays non-radiatively to the lower branch within the ra- 
diative lifetime of ^E, and that there exists no commensurate non-radiative 
process in the opposite direction from the lower branch to the upper branch, 
thereby making the emission from the lower branch much more intense than 
from the upper branch. The dominance of the non-radiative decay from the 
upper branch to the lower branch is due to the presence of the large stress 
splitting and the low temperature conditions of Davies' and Hamer's mea- 
surement. Consequently, it appears that in using such high stress and low 
temperature, Davies and Hamer observed just the small Jahn- Teller effect 
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Figure 24: Polarised luminescence spectra of NV~ as taken by Davies and 
Hamer [7T] at 17 K using a high density ensemble and an applied stress of 1.4 
GPa along the [001] direction. The left panel shows that the ZPL has been 
split into upper (line a) and lower (line b) branches, which are expected to be 
purely vr-polarised (parallel to stress axis) and cr-polarised (perpendicular to 
stress axis), respectively. The right panel shows that the polarisation of the 
lower branch (solid curves) is a mixture of both polarisations, which Davies 
and Hamer explained by the presence of random crystal strain. The left 
panel represents the polarisation of the luminescence band once corrected for 
the depolarising effects of random crystal strain inferred from the lower ZPL 
branch. 



that arises from the time-independent mixing of the vibronic levels of 
due to relatively strong interactions with pseudo-local phonon modes of E 
symmetry. This is unlike the recent experiments performed by Fu et al [SB] 
and Kaiser et al |142] using single centres at low temperatures (similar to 
those employed by Davies and Hamer) with much smaller strain splittings 
(< 100 GHz), where rapid two-way non-radiative transitions between the 
ZPL branches occurred. As will be now discussed, Fu et al and Kaiser et al 
observed the Jahn- Teller effect that arises from both the time-independent 
vibronic coupling due to E symmetric pseudo-local modes and the rapid 
two-way transitions between the ZPL branches that are mediated by the 
continuum of non-local E symmetric phonon modes. 

In their low temperature (~ 4 K) single centre study. Kaiser et al |142] 
performed two measurements to establish the polarisation selection rules for 
the absorption and emission of the NV~ optical ZPL. The centres they stud- 
ied exhibited small strain splittings (3-10 GHz). Their first experiment con- 



56 



sisted of the polarised resonant excitation of one of the ZPL branches and 
the detection of the fluorescence intensity as a function of excitation polar- 
isation. Kaiser et al observed that the ZPL branches were orthogonal in 
excitation polarisation with ~100% contrast, implying that the ZPL absorp- 
tion selection rules were almost perfect. Their second experiment consisted 
of the resonant excitation of one of the ZPL branches and the detection of the 
fluorescence polarisation. Kaiser et al observed that the emitted polarisation 
contrast was < 60% for all of the centres they investigated and suggested 
that the variance of the emission contrast between centres is likely due to 
variations in the local environment between centres. When this second result 
is combined with the observation of Fu et al [88] (in their similar low temper- 
ature single centre study) that the emission polarisation of each ZPL branch 
is well polarised with ~ 100% contrast, they imply that the ZPL emission 
selection rules are perfect and that the emission depolarisation is occurring 
due to non-radiative population transfer between the ZPL branches within 
the radiative lifetime of ^E. This conclusion is consistent with the preceding 
discussion regarding the luminescence spectra of Davies and Hamer [7T] as 
well as the observation of orbital averaging of the fine structure at room 
temperature [87] discussed in the previous section. 

Fu et al also measured the temperature dependence of the emission po- 
larisation contrast of the ZPL branches of single centres with various strain 
splittings (8-81 GHz) under polarised off-resonance (532 nm) excitation in 
the temperature range 10-50 K. In agreement with the resonant excitation - 
fluorescence detection measurements of Kaiser et al, Fu et al observed < 60% 
polarisation contrast at all temperatures. Note the subtle, but complemen- 
tary, difference between the two investigations is that Fu et al excited into 
the absorption band and measured just the emission polarisation of the ZPL 
branches, whereas Kaiser et al excited into the ZPL and measured the polar- 
isation of the entire emission band. Furthermore, Fu et al observed that the 
emission polarisation contrast depended significantly on temperature, such 
that by 40 K, the contrast had reduced to ~ 0% and was independent of the 
excitation polarisation (refer to figure [5]). Fu et al noted that the temperature 
dependence of the emission polarisation contrast differed between individual 
centres and demonstrated a correlation with strain. In order to explain their 
observations, both Kaiser et al and Fu et al used a three-level rate equation 
model to describe the dynamics of the optical transitions between and ^E 
and the non-radiative transitions between the fine structure orbital branches 
of ^E. Fu et al provided a more general model and that model is depicted in 
figure [231 The key difference in the models of Fu et al and Kaiser et al was the 
inclusion of the excitation cross transition with rate 7a by Fu et al. Kaiser et 
al did not include this cross transition because it represents imperfect ZPL 



57 



absorption selection rules for their case of resonant excitation, which is con- 
tradictory to their first observation. Fu et al included the cross transition to 
account for the possibility of imperfect absorption selection rules when the 
optical transition is excited off-resonance. For individual centres with strain 
splittings of ~ 8 GHz, Fu et al found that the non-radiative rate between the 

orbital branches depended on temperature as the fifth power of temper- 
ature b/c oc T^, and that the cross transition ratio was 7 ~ 40%. The tem- 
perature dependence of the non-radiative rate is indicative of a two-phonon 
Raman transition between the orbital branches involving phonons with fre- 
quencies much greater than the strain splitting of the branches. Given the 
dependence of the non-radiative transitions between the branches, the cross 
transition rate for off-resonance excitation governs the temperature depen- 
dence of the polarisation contrast at very low temperatures, such that there 
remains < 60% polarisation contrast in the limit of zero temperature. Fu 
et al attributed the origins of the cross transition and their consequence of 
imperfect polarisation contrast at very low temperature to the Jahn- Teller 
effect, but did not provide a thorough explanation. 

Based upon the observations of Davies and Hamer [71], up to ~ 10% of 
the > 40% depolarisation observed by Fu et al at very low temperature can 
be attributed to temperature independent cross transitions allowed by the 
time- independent Jahn- Teller effect arising from the couplings of the vibronic 
levels of that are induced by E symmetric pseudo-local modes. A small 
amount of the remaining unexplained depolarisation at low temperature can 
be attributed to direct one-phonon transitions between the ZPL branches 
that was neglected in the rate equation model of Fu et al. However, it is 
evident that a significant portion of low temperature depolarisation remains 
unexplained and is an important issue in the understanding of the NV~ centre 
to resolve. One possible avenue of explanation is a strain dependence of the 
time-independent Jahn- Teller effect, such that in low strain conditions, the 
depolarisation arising from the time-independent Jahn- Teller effect is much 
greater than the ~ 10% depolarisation observed by Davies and Hamer in 
high strain conditions. 

The rapid non-radiative transitions between the ZPL branches also con- 
tribute to the homogeneous broadening of the ZPL. In their low temperature 
single centre study, Fu et al [88] performed the first direct measurements of 
the homogeneous ZPL width as a function of temperature. Fu et al's sin- 
gle centre measurements were unlike the ensemble measurements performed 
by Davies [171j and later by Hizhnyakov et al [141j . who only obtained the 
homogeneous ZPL width indirectly via a deconvolution with the ensemble 
inhomogeneous lineshape. Fu et al found that the upper ZPL branch ho- 
mogeneously broadens also as the fifth power of temperature, which corre- 



58 




7a 



'A2^ 

Figure 25: (a) The Jahn- Teller model of fluorescence depolarisation and ho- 
mogeneous ZPL broadening of the NV~ optical transition provided by Fu et 
al [SB]. The three electronic levels denote and the two strain split low 
temperature fine structure branches of that are orthogonally polarised. 
The solid arrows denote optical absorption and emission transitions with 
rates a and c, respectively. The optical excitation is taken to be ^/-polarised, 

such that the dashed ( ) arrow denotes the potential cross transition that 

populates the x-polarised branch due to the modification of off-resonance 
excitation selection rules by the Jahn- Teller effect. The chain (— ■ — ) arrows 
denote the non-symmetric phonon transitions between the orbital branches 
with decay rate b and a Boltzmann prefactor e"''/^-^^ excitation, where 6± 
is the strain splitting between the two branches. 



lates with the two-phonon Raman transitions between the orbital branches 
of that govern the depolarisation of the optical emission. Given that 
the symmetric mode model of Davies |171j yields a significantly different 
homogeneous broadening of the ZPL, it may be concluded that the fit that 
Davies obtained via deconvolution with a temperature dependent inhomoge- 
neous lineshape was fortuitous and that the symmetric mode model cannot 
explain the homogeneous broadening of the ZPL. The observation of a 
homogeneous broadening is supported by the ultrafast photon echo measure- 
ments of ensembles at low temperature performed by Lenef et al [156] , who 
observed a dephasing of the upper ZPL branch. However, Lenef et al 
observed a dephasing of the lower ZPL branch (not measured by Fu et 
al) , which is supported by the indirect ensemble measurements of the inho- 
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mogeneously broadened ZPL performed by Hizhnyakov et al |141] . Lenef 
et al speculated that the dephasing arose from the combination of the 
Jahn- Teller effect and large crystal strain, but did not provide any specific 
details. Alternatively, Hizhnyakov et al speculated that the broadening 
is due to a combination of a strong linear and weak quadratic Jahn- Teller 
effects, but there is no other evidence of a strong linear effect. Consequently, 
the homogeneous broadening of the lower ZPL branch still requires further 
experimental investigation, preferably by the direct means employed by Fu 
et al. 

In spite of the conclusive empirical evidence of the presence of the Jahn- 
Teller effect in ^E, there has not yet been a thorough attempt at a Jahn- Teller 
explanation of the differences in the NV~ optical absorption and lumines- 
cence bands. In their ah initio calculation, Gali et al [521 identified Ai and 
E pseudo-local modes with vibrational energies 76.9 meV and 63.9 meV as- 
sociated with "^E, respectively, and suggested that if the Jahn- Teller effect 
was present in ^E, then optical transitions would be allowed from the ground 
vibronic level of '^A2 to the first vibronic levels of ^E corresponding to these 
modes (refer to figure [2^ . thereby yielding a double one-phonon peak in 
absorption. Whilst this conclusion was correct, Gali et al failed to identify 
that if this was the case then the modified selection rules would induce a 
similar doubling of the one-phonon peak in emission arising from the first 
vibronic levels of corresponding to the same Ai and E modes, which 
have ground state vibrational energies 77.0 meV and 63.3 meV, respectively. 
Gali et al's assertion is thus clearly incorrect. Another ah initio calculation 
performed by Zhang et al [70] provided additional evidence for a weak Jahn- 
Teller effect in ^E by finding that the ^E electronic energy is reduced by ~6.4 
meV, if the symmetry of the centre is relaxed from C^^, to Cih- Zhang et al 
identified that the mode with a displacement coordinate equivalent to the 
symmetry lowering relaxation had an energy much greater than ~6.4 meV. 
Consequently, this result implies a weak dynamic Jahn- Teller effect where 
the zero-point energy is greater than the Jahn- Teller distortion energy |70] . 
In summary, the empirical and ah initio evidence indicates the presence of a 
weak dynamic Jahn- Teller in ^E and the thorough theoretical treatment of 
such an effect remains an outstanding issue of the NV~ centre that needs to 
be addressed. 

4.5 The NV infrared vibronic band 

The infrared vibronic band of NV~ differs to the optical band in several 
ways. Since the infrared band occurs between the intermediate singlet states 
{}Ai and ^E), it has so far only been observed in fluorescence via the ex- 
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Figure 26: The Jahn- Teller model of the NV~ optical transition proposed by 
Gali et al [69]. The vibronic levels are denoted by the electronic state (left) 
and vibrational occupation (right) of the Ai and E modes that Gali et al cal- 
culated to have vibrational energies 77.0 meV and 63.3 meV associated with 
^A2 and energies 76.9 meV and 63.9 meV associated with ^i?, respectively. 
The solid arrows indicate optical transitions originating from the ground vi- 
bronic levels of each electronic state that are allowed in the absence of the 
Jahn- Teller effect and the dashed arrows indicate the transitions that are 
allowed as a consequence of the effect. The transitions will result in a double 
one-phonon peak in both the absorption and luminescence bands in conflict 
with experiment [7T]. Gali et al did not identify the dashed transition in 
emission. 



citation of the centre's optical transition and the non-radiative transfer of 
population to the upper singlet ^Ai [771 1165[ I164j . Consequently, the emis- 
sion band has only been measured to date and there is no observation of 
the absorption band available for comparison. The emission band (refer to 
figure [TS]) exhibits a one-phonon peak with an energy of 42.6 meV from the 
ZPL and a temperature independent total Huang- Rhys factor of S" ~ 1 [77] . 
The infrared emission band is thus very different to the optical band, where 
several phonon peaks occur at energies that are almost perfect multiples of 
the one phonon energy (~ 67 meV) [7T]. A further difference between the 
infrared and optical bands is the temperature independence of the infrared 
band [77 ] 1164] . Unlike the optical band, the features of the infrared emission 
band have been observed to be approximately temperature independent [77] 
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and the infrared ZPL has been observed to not significantly broaden or shift 
in energy with temperature |164] . Applying the symmetric mode model, the 
temperature independence of the infrared emission band together with its 
small Huang-Rhys factor imply that interactions with symmetric modes do 
not play a significant part in forming the infrared band. This conclusion is 
consistent with a simple molecular model, where since ^Ai and belong to 
the same MO configuration afe"^, there is no axial shift of electronic charge 
involved in the infrared transition and therefore, no excitation of symmetric 
modes. Hence, it appears that the infrared band predominately arises due 
to interactions with non-symmetric modes. 

As discussed in the previous section and consistent with the preceding 
conclusion, the Jahn- Teller effect has been observed in via a vibronic 
LAC under uniaxial stress by Manson et al [165] . The presence of the LAC 
can be explained by first identifying that in the absence of stress, the Jahn- 
Teller interaction of ^E with a single pair of E symmetric vibrational modes 
gives rise to vibronic levels of E, Ai and A2 symmetry, as depicted in figure [271 
|165] . The ground vibronic state of E symmetry splits in the presence of non- 
symmetric strain and interacts quadratically with the next vibronic level of 
Al symmetry |165j . At large enough stress, the quadratic interaction between 
the vibronic levels will result in a LAC. Importantly, the Ai vibronic level is 
not observed in the infrared spectra until the LAC is reached, since transitions 
from ^Ai to the Ai and A2 vibronic levels of ^E are forbidden unless there 
is significant stress mixing of those vibronic levels with the ground vibronic 
level of ^E. This behaviour enabled the identification of the Ai vibronic level 
and furthermore, implied that the lower vibronic levels of^E are energetically 
well spaced, as there were no additional LACs observed [165] . 

From their stress measurements, Manson et al inferred that the Ai vi- 
bronic level is ~ 14.3 meV above the ground vibronic level of ^E at zero 
stress. The energy of the E vibronic level above the ground vibronic level 
was directly identified from the emission band as being at 42.6 meV, since 
the Jahn- Teller coupling of the lower E vibronic levels of ^E allows radia- 
tive transitions between them and ^Ai with and without the application of 
stress [165] . Consequently, the A2 vibronic level must exist in the range 
14.3 — 42.5 meV with an energy that is sufficiently greater than 14.3 meV 
in order to not be observed in the uniaxial stress measurements. As con- 
cluded by Manson et al, the observed vibronic levels of ^E can only occur 
if both linear and quadratic Jahn- Teller interactions are significant. In ad- 
dition, since a vibronic LAC could only occur in ^E (as ^Ai cannot interact 
with non-symmetric modes) and that it must be on the high energy side of 
the ^E ground vibronic level (as there are obviously no vibronic levels below 
the ground vibronic level) , the observation of the vibronic LAC enabled the 
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Figure 27: (a) The vibronic structure related to the infrared band of NV~ 
|165] . The ground (0) vibronic level of ^Ai and the ground and first excited 
(1^;) vibronic levels of corresponding to the occupation of a pair of E 
symmetric modes in the absence of the Jahn-Teller effect is depicted on the 
left hand side. The vibronic splittings due to linear and quadratic Jahn- 
Teller interactions are depicted in the centre and the interactions of the split 
vibronic levels in the presence of non-symmetric stress is depicted on the 
right hand side. Solid arrows originating from ^Ai denote allowed infrared 
transitions in the absence of the Jahn-Teller effect (left) and in the presence of 
the effect (right). The dashed arrow denotes the additional infrared transition 
to the Ai vibronic level that is allowed in the presence of stress. The solid 
arrows in the centre denote the known vibronic energies |165j . 



unambiguous ordering of the intermediate ^Ai and ^E singlets (as the LAC 
is observed in the IR emission band) |165] . 

4.6 The NV^ optical vibronic band 

The optical band of NV'' has aspects that are similar to both the optical and 
infrared bands of NV~. Since the NV° '^E{a\e)—'^Ai{aie^) optical transition 
involves a similar change in MO configuration to the NV~ optical transition, 
interactions with symmetric modes contribute significantly to the NV° optical 
band. The low temperature total Huang-Rhys factor of the NV° optical 
band is S* ~ 3.3 and the ZPL energy, intensity and width are dependent 
on temperature |176j . Additionally, the Jahn-Teller effect is observed to be 
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present in "^E and it produces an analogous vibronic LAC as observed in 
of NV~ under uniaxial stress [70]. However, the energy ordering of the Ai 
and A2 vibronic levels of '^E is unknown since the orbital symmetry of the 
optically excited of NV° has not been established. The lowest energy A 
vibronic level of "^E occurs at ~ 12 meV, which is a very similar energy to 
the Ai level of ^E [TB]. Although there does not currently exist a detailed 
analysis of the NV~ infrared band, from first impressions it is expected that 
there will be significant similarity between the Jahn- Teller interactions in '^E 
of NV° and in of NV". 

4.7 The temperature dependence of the fine structure 
of the NV~ level 

Apart from their presence in the optical and infrared bands of NV~ and 
NV°, potential vibronic effects have also been observed in the temperature 
dependence of the NV~ ground state fine structure. In reference to the 
NV~ ground state spin-Hamiltonian and potential for interactions with strain 
fields, Acosta et al [29], I177j observed that the fine structure splitting Dgs 
between the = and = ±1 spi n sub-l evels varied with temperature 
and that the strain parameter 5± = + Sy that splits the = ±1 spin 
sub-levels in the presence of transverse strain did not vary with temperature. 
Acosta et al attempted to describe the temperature dependence of Dgg by a 
model based purely on the thermal expansion of the crystal lattice, but were 
only partially successful. They ultimately concluded that additional vibronic 
effects beyond lattice expansion must be responsible for the temperature 
dependence of Dgs. The identification and description of these vibronic effects 
remain outstanding issues concerning NV~. 

By comparing the temperature dependence of Dgs and that of the NV~ 
optical ZPL energy it is evident that their dependence is very similar, which 
suggests that quadratic symmetric mode coupling could be responsible for the 
temperature dependent shift of Dgs, just as it was for the optical ZPL. This 
proposal is supported by the observation that the strain splitting 6± does 
not change with temperature, thereby implying that interactions with non- 
symmetric modes are not sufficient to reduce the splitting as they do to the 
analogous splitting of ^E in the process of orbital averaging. Consequently, 
it is clear that there is potentially a great deal of information concerning 
vibronic interactions in the NV centre available through the study of the 
NV~ ground state temperature dependence. 

This section has been principally focussed on the vibronic interactions 
that occur within a given electronic state and give rise to directly observable 
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vibronic structure. In the next section, the key observations of the vibronic 
interactions between electronic states that give rise to non-radiative transi- 
tions will be discussed. As will be shown, these transitions have a significant 
role in the optical and spin dynamics of NV~ and are central to understand- 
ing the properties of optical spin-polarisation and long lived spin coherence. 

5 Optical and spin dynamics 

Manson et al [IZl [52] were the first to distinguish the fast and slow compo- 
nents of the NV~ optical dynamics through dual optical pulse measurements. 
The fast dynamics were attributed to intrinsic radiative and non-radiative 
transitions between the electronic states of NV~ and the slow dynamics were 
attributed to extrinsic transitions involving photoconversion to NV". The 
current understanding of the intrinsic and extrinsic aspects of the optical 
dynamics of NV~ are depicted in figure [281 The intrinsic dynamics are re- 
sponsible for optical spin-polarisation and readout. The extrinsic dynam- 
ics are responsible for the establishment of a dynamical equilibrium of the 
charge states when both are simultaneously excited optically as well as the 
photoconversion of NV~ into NV° when just NV~ is excited. Due to their 
significantly different timescales, the extrinsic and intrinsic optical dynamics 
of NV~ can be treated individually in most conditions. 

5.1 The extrinsic optical dynamics of NV~ 

It is currently understood that the extrinsic dynamics involve both radiative 
and non-radiative transitions. The radiative transitions between the charge 
states are believed to occur between their optically excited levels via interme- 
diate valence/conduction band derived levels |lU7j . The radiative transitions 
between the charge states therefore correspond to photoconversion processes 
that are quadratic in the optical excitation intensity. Furthermore, since 
the radiative transitions do not explicitly involve electron donors/acceptors, 
but simply the inherent conduct ion/ valence bands of the crystal, the ra- 
diative transitions are expected to be independent of the concentrations of 
donors/acceptors in the crystal. Differing from the radiative transitions, the 
non-radiative transitions occur directly between the electronic states of the 
charge states. Consequently, the non-radiative transitions must correspond 
to the tunneling of an electron between the NV centre and proximal electron 
donors/acceptors. Such tunneling will be highly dependent on the distance 
between the NV centre and the donors/acceptors, which will be itself statis- 
tically dependent on the concentrations of donors/acceptors. Apart from the 
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Figure 28: The current understanding of the extrinsic and intrinsic aspects 
of the optical dynamics of NV^. Sohd arrows denote radiative transitions, 

dashed ( ) arrows denote non-radiative transitions within a charge state 

and chain (— ■ — ) arrows denote non-radiative transitions between charge 
states. Electronic states of NV° are depicted on the left hand side, whereas 
electronic states of NV~ are depicted on the right hand side. The optical ZPL 
energies of each charge state indicated. The upper most and central 

electronic state represents the continuum of states derived from the occu- 
pation of the conduction band. The conduction band states are believed to 
form an intermediate step in the quadratic photoconversion process between 
the NV charge states. In addition to the non-radiative transitions between 
the ground states of the two charge states, only non-radiative transitions 
that originate from their optically excited states are also depicted. Note that 
the final states of these transitions are unknown (as indicated by question 
marks) and that other non-radiative transitions between the charge states 
can potentially occur between any of the electronic states. 
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non-radiative decay from the ground electronic state of the unstable charge 
state (NV° in the case of figure [25]) to the ground electronic state of the 
stable charge state (NV~ in the case of figure [28]) in the absence of optical 
illumination |120] . the other non-radiative transitions originate from the op- 
tically excited levels of each charge state. Thus, the non-radiative transitions 
correspond to photoconversion processes that are linear in optical excitation 
intensity [S3]. 

The first indication that quadratic photoconversion processes occur in the 
NV centre was identified by Redman et al |178j . who observed a weak pho- 
toconductive current consistent with photoionisation occurring from a satu- 
rated excited electronic state of NV~. This observation was later supported 
by the single centre studies [1621 1179] that attempted to fit the observed NV~ 
fluorescence autocorrelation function for different excitation intensities using 
the three level rate equation models that were discussed in section [3] In order 
to correctly fit the autocorrelation function, the transition rates between 
and the long lived metastable state of the model had to be considered to be 
proportional to the optical excitation intensity. It is now understood that 
under the strong excitation intensities used in those studies, the long lived 
metastable state of the three level models effectively represented NV° and 
thus, the intensity dependent transitions to and from the metastable state 
were in fact representative of the radiative transitions between the charge 
states [S2]. However, the conclusive evidence that quadratic photoconver- 
sion processes occur in the NV centre was provided by Waldherr et al )1U7] . 
who observed that when NV~ was continuously excited on resonance (1.945 
eV) the fluorescence would decay exponentially and if, when the emission 
extinguished, an off-resonance (2.32 eV) excitation was instead applied, the 
NV~ fluorescence would recover exponentially. Waldherr et al found that 
the decay and recovery rates were quadratic in excitation intensity for low 
intensities and linear at higher intensities. Furthermore, they identified by 
observing the changes in NMR spectra that under resonant /off- resonance op- 
tical excitation, the centre was being converted from NV~ to NV° and back 
again. Hence, this identification and the dependence of the decay/recovery 
rates on excitation intensity conclusively demonstrated that radiative tran- 
sitions between the charge states occur in each direction from the optical 
excited states of NV~ and NV°, such that the excited states become satu- 
rated at sufficiently high excitation intensities )107] . 

Direct observations of the linear photoconversion processes were made 
by Manson et al [M], who found that the ratio of the ZPL emission in- 
tensities of NV~ and NV'' follows a saturation curve when excited by an 
optical energy (2.32 eV) that is within the optical absorption bands of each 
charge state. The implication of the occurrence of linear processes is that 
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the conversion between the charge states must proceed in both directions via 
a non-radiative decay from the optically excited states of the charge states. 
As discussed, such non-radiative decay requires the tunneling of an electron 
between the NV centre and proximal electron donors/acceptors, whose dis- 
tances from the NV centre are statistically determined by the concentrations 
of donors/acceptors [53]. Consistent with this conclusion, the ensemble mea- 
surements of Manson et al were performed using type lb diamond samples, 
which contain high substitutional nitrogen (electron donor) concentrations. 
The conclusion also implies that linear processes are expected to occur in NV 
centres in type lb nanodiamonds, where there are significant concentrations 
of both substitutional nitrogen and surface traps, however no studies have 
been conducted to date to confirm. 

Although the extrinsic dynamics appear to be reasonably established, the 
understanding of their origins is certainly not. To date, there has not been 
a thorough microscopic model developed to explain the various radiative 
and non-radiative transitions between the charge states beyond simple spec- 
ulations of the role of conduction and/or valence band states and electron 
donors/acceptors. There has also been no attempt to relate the dynamics 
to the concentrations of the charge states, spectral stability of the ZPLs or 
material properties. Consequently, the microscopic origins of the extrinsic 
dynamics and the stable charge state of the centre remain outstanding issues 
that need to be addressed. 

5.2 The intrinsic optical dynamics of NV~ 

Turning now to the intrinsic dynamics of NV~, the intrinsic dynamics are 
most efficiently discussed using a six level model, such as the one depicted 
in figure 1211 The six level model is an accurate representation of the NV~ 
electronic structure at room temperature, however more complicated models 
are required when magnetic fields are applied such that the spin sub-levels 
interact or when operating at low temperatures with resonant excitation. 
Measurements of the transition rates of the six level model are contained 
in table [TDi Note that transitions between the spin sub-levels of each triplet 
state are not included in the six level model of figure [21] since they are several 
orders of magnitude slower than the transitions of the optical dynamics. 
The transitions between the spin sub-levels of each triplet are considered 
to be aspects of the NV~ spin dynamics, which will be discussed shortly. 
The other notable absences from the six level model of figure [21] are the 
transitions directly from to and from ^Ai to ^A2. These transitions 
are of sufficiently high energy (> 1.2 eV) that radiative decay will be expected 
to dominate any many-phonon non-radiative decay. Since there has been no 
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detection of these transitions in the fluorescence spectra of NV~ to date, they 
are beheved to be neghgible compared to the other transitions of the intrinsic 
optical dynamics. 
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Figure 29: The six level model of the intrinsic dynamics of NV~. Electronic 
levels are labelled by their level index (1-6) on the left hand side and spin 
sub- level on the right hand side (triplet states only). The solid arrows denote 

transitions that conserve spin quantum numbers, dashed ( ) arrows denote 

cross transitions between triplet spin sub- levels and chain (— ■ — ) arrows 
denote intersystem crossing transitions between triplet and singlet states. 



Collins et al [180] were the first to measure the lifetime of ^E and ob- 
tained 12.9±0.1 ns for NV~ ensembles in natural type lb diamond. Collins 
et al stated that the lifetime was independent of temperature in the range 
77-700 K, which has been subsequently supported by another study [164] . 
but seemingly contradicted by a recent third study performed by Toyli et al 
|30j . The inconsistency of Toyli et al's lifetime measurements has not been 
resolved and it is an important issue to be pursued. For the purposes of this 
discussion, the lifetime will be considered to be temperature independent. 
The first single centre studies [H [2] reported the observation of biexponential 
decay of the '^E fluorescence with time instead of the monoexponential decays 
observed by Collins et al and other previous ensemble measurements [156] . 
Using three level models, the single centre studies attributed the biexponen- 
tial decay to the repopulation of ^E from the long lived metastable state 
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Table 10: Experimental measurements of the NV~ excited state lifetimes tab- 
ulated by reference. The lifetimes are defined in terms of the transition rates 
kij of the six level model depicted in figure by = 1/ J2j hj- Uncertainty 
ranges are included if provided by sources. 



Reference 


T(K) 


T-i (ns) T4 (ns) 


T5 (ns) 


Te (ns) 


Colhns^ |180j 


77-700 




12.9 ±0.1" 




Lenef^ [T^ 


77 




12.96±0.14" 




Manson^ [17] 


Room 


300'^ 






Batalov* [96] 


Room 




12.0 


7.8 


Acosta" [29] 


4.4-70 


462±10^ 0.9±0.5 


10.0±0.3" 




Robledo" [97] 


8 




10.9 




Robledo^ [181] 


300 


371±6^ 


13.7±0.1 


7.3±0.1 


Beveratos"'" [TS2] 


Room 




25" 




Beveratos"''" [37] 


Room 




23" 




Neumann*'" [80] 


Room 




23 


12.7 



"did not distinguish triplet spin sub- levels/ singlets; ''low T<5 K lifetime 
'^ensemble; '^single centre; "nanodiamond 



[H [2]. This is now known to be incorrect and the existence of a biexponen- 
tial decay is the result of the spin sub-levels of having different lifetimes 
due to differing non-radiative decay rates to the intermediate singlet states. 
Batalov et al [96] were the first to measure the individual lifetimes of the 
spin sub-levels in bulk diamond. The lifetimes obtained by Batalov et al are 
largely consistent with the later measurements performed by Robledo et al 
|181j . The ^E lifetimes approximately double in nanodiamond crystals that 
are much smaller than the fluorescence wavelength. This change is simply 
due to the reduction of the radiative emission rate induced by the decrease of 
the effective refractive index of the medium surrounding the NV centre. Bev- 
eratos et al |182j were the first to measure the ^E lifetime of a single centre in 
such a nanodiamond and observed a monoexponential decay corresponding 
to a lifetime of 25 ns. Neumann et al subsequently measured the individual 
lifetimes of the ^E spin sub-levels for a single centre in a nanodiamond. 

The equality of the ^E lifetime obtained from observations of monoexpo- 
nential fluorescence decay and the = spin sub-level lifetime obtained 
from observations of biexponential decay, strongly suggests that the experi- 
mental conditions in which monoexponential decay was observed were such 
that the centre had been spin-polarised into the = spin-projection prior 
to the recording of the fluorescence decay curve [17] . Such preliminary spin- 
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polarisation could have occurred if optical pulses that were too intense or 
long or insufficiently spaced were used to conduct the lifetime measurement. 
Furthermore, the temperature independence of the lifetime of the rus = 
sub-level, even up to temperatures (700 K) where the highest energy phonons 
(~160 meV) of the diamond lattice are appreciably occupied, unambiguously 
implies that the non-radiative decay from the = sub-level is negligible 
compared to the radiative decay. This conclusion is supported by the low 
temperature optical Rabi experiments of Batalov et al [96] , in which nearly 
pure radiative relaxation and dephasing were observed. In the absence of 
any strong spin-dependent interactions between other electronic states and 
and ^A2, which have not been detected to date in any of the uniaxial 
stress, magnetic and electric field studies, the radiative decay rates of will 
be spin independent |17]. Thus, the difference between the lifetimes of the 
^E rris = and rUs = ±1 spin sub-levels is due to the m<j = ±1 sub-levels 
undergoing greater non-radiative decay, which ultimately results in the spin- 
dependent fluorescence of NV~ [17] . Using the observed sub- level lifetimes, a 
simple expression can be derived that relates the non-radiative contributions 
to the decay of each spin sub-level 

n 1 + 70 (^Q) 



Ts 1 + 7± 



where is the ratio of the non-radiative decay rate of sub-level to the 
^E radiative decay rate. In the limit 70 ~ 0, the lifetime measurements of 
Batalov et al result in the estimation 7-1-1 f» 0.54, which increases if 70 > 0. 
Hence, it may be concluded that the non-radiative decay rate out of the ^E 
nis = ±1 spin sub-levels is at least 50% of the ^E radiative decay rate [17] . 
This conclusion is consistent with the lifetime measurements of single cen- 
tres in nanodiamond crystals performed by Neumann et al. Such a significant 
contribution of non-radiative transitions to the decay of the ^E rris = ±1 spin 
sub-levels should ensure that their lifetime displays an observable tempera- 
ture dependence. However, in the recent lifetime measurements of Toyli et 
al [30], the lifetimes of the m<, = ±1 spin sub- levels were observed to exhibit 
a small temperature reduction and the rris = spin sub-levels were observed 
to undergo a much larger temperature reduction. These observations are 
inconsistent with the preceding discussion and pose an issue that needs to 
be resolved. 

At low temperatures, each of the ^E rrig = ±1 fine structure levels are 
expected to have different lifetimes because of their different spin-orbit sym- 
metry (refer to figure [T7I) . Simple symmetry arguments imply that the £'(±1) 
fine structure level can only decay to ^Ai through the emission of E sym- 
metric phonons, whereas the ^i(±l) level can decay directly to ^Ai through 
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the emission of Ai phonons or in a two step process via the £'(±1) level 
through the emission of E phonons. The y42(±l) fine structure level can 
only decay to ^Ai in a two step process via the -£^(±1) level through the 
emission of E phonons. Consequently, based purely on the possible decay 
paths for each rris = ±1 fine structure level, /li(±l) is expected to have 
the shortest lifetime, /l2(±l) the longest lifetime and £'(±1) an intermediate 
lifetime. As mentioned in section [31 such simple lifetime arguments were 
used to determine the energy ordering of the Ai and A2 fine structure levels 
through the inspection of the low temperature excitation spectra of the ^E 
fine structure (refer to figure 1301) [TO] . At room temperature, orbital averag- 
ing is expected to lead to the single lifetime of the = ±1 spin sub-level in 
the six level model (refer to figure [2^ . A detailed model of the non-radiative 
decay mechanisms of the low temperature fine structure levels and their av- 
eraging at room temperature has not yet been developed and remains an 
outstanding issue concerning NV~. 

The radiative non spin-conserving transitions between spin sub-levels of 
the triplet states with different spin-projections occur due to the small mixing 
of the spin-projections in ^E arising from spin-spin interaction (parameter 
X-^g of the spin-Hamiltonian ([3])). Given the observed magnitude of the spin- 
spin interaction (refer to table [3]) being an order of magnitude smaller than 
the low temperature low strain fine structure splittings of ^E, the rates of 
the cross transitions are expected to be less than 1% of the spin-conserving 
radiative transition rate. This expectation has been supported by recent low 
temperature measurements using resonant excitation of the fine structure 
levels [To]. As discussed in section [3l in the presence of higher strain, the 
spin interaction causes the significant spin mixing in the lower branch of 
the ^E low temperature fine structure that gives rise to optical spin-flip 
transitions. Therefore, the relative strength of the cross transitions are a 
function of the strain present at the NV centre. No model has yet been 
developed to describe the strain dependence of the cross transitions at room 
temperature. Regardless, it is expected that the radiative cross transitions 
play a minor role in the dynamics that lead to optical spin-polarisation, since 
the transitions are equal and thus, do not result in the preferential population 
of one particular spin-projection [171 1165] . 

As discussed in section [3l the lifetimes of the intermediate singlet states as 
well as their temperature dependence have been measured using NV~ ensem- 
bles by Acosta et al [29]. The lifetime of ^Ai was observed to be temperature 
independent in the temperature range 4.4—70 K. Since it is currently believed 
that the non-radiative decay between the singlet states dominates the weak 
infrared decay [77] , the observed temperature independence of the lifetime of 
^Ai implies that the non-radiative decay out of ^^41 occurs via high energy 
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Figure 30: Low temperature optical excitation spectra of the NV~ fine 
structure obtained from a single centre using continuous microwave driving 
of the ground state magnetic resonance [10]. The fine structure levels corre- 
sponding to the spectral lines are denoted by their spin-orbit symmetry and 
spin-projection. The splitting of the central -E'(O) level into two lines signifies 
that there is a small strain present in the centre of ~ 3 GHz. The lower 
energy E{0) line is slightly less intense than the upper E{0) line due to the 
small enhancement of the spin-spin coupling of the lower energy E{0) branch 
with the E{±1) level in the presence of strain. The relative intensities of the 
upper two lines were used to infer their spin-orbit symmetry. The least in- 
tense line corresponds to the Ai(±l) level since it has the most non-radiative 
decay paths. The highest energy line is the most intense out of the rris = ±1 
lines and thus, corresponds to the A2(±l) level since it has the least number 
of non-radiative decay paths. 



phonons, whose occupations do not appreciably change in the temperature 
range 4.4 — 70 K. For example, the thermal occupation of a phonon with 
energy 67 meV is < 3 x 10~^ at 70 K. The lifetime of ^E, which is expected 
to undergo pure non-radiative decay to the ground state, on the other 
hand displays a significant temperature dependence [221 11^1] • This signifies 
that the non-radiative decay between the singlet states occurs via different 
electron-phonon interactions than the decay from ^E to ^A2. Acosta et al 
found that the temperature dependence of the ^E lifetime could be modeled 
by the emission of a phonon of energy 15(1) meV. The later single centre 
lifetime measurements of Robledo et al [181] were consistent with those of 
Acosta et al and obtained a phonon energy of 16.6 ± 0.9 meV. 

Robledo et al concluded that the phonon energy was indicative of a third 
electronic state ~ 16 meV below the lower energy electronic state involved 
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in the infrared transition. Motivated by the ab initio calculations of Ma et 
al [HB], Robledo et al speculated that there existed three intermediate sin- 
glet states and that the lower energy electronic state involved in the infrared 
transition was ^Ai, which decayed via a ~ 16 meV phonon to an even lower 
energy ^E. Given that Ma et al's calculation has been shown to be inconsis- 
tent with the observations of Manson et al |165] and that the ordering of the 
singlets involved in the infrared transition has been established as depicted 
in figure [291 the speculation of Robledo et al is also brought into question. 
Given the detection by Manson et al of an Ai vibronic level ~ 14.3 meV 
above the ground vibronic level of ^E, the equivalence of the energy of the 
vibronic level and the observed phonon energies strongly suggests that the 
decay from ^E involves the Ai vibronic level instead of an unknown elec- 
tronic state. Using this assertion, the non-radiative decay from ^E to ^A2 
most likely originates from both the ^E ground and Ai (occupied in thermal 
equilibrium) vibronic levels and involves the emission of multiple high energy 
phonons that have no observable temperature dependence in the tempera- 
ture ranges investigated by either Acosta et al or Robledo et al |165j . This 
explanation is consistent with the vast majority of theoretical models and 
empirical observations to date, which do not indicate the presence of a third 
intermediate singlet state. However, this explanation is yet to be developed 
in any detail and must include the consideration of the Jahn- Teller effect 
in ^E as well as the spin dependent electronic interactions that enable the 
non-radiative decay to ^^2. 

Using a sequence of ps optical pulses combined with the microwave control 
of the ground state electronic spin, Robledo et al also measured the ratio of 
the decay rates from ^E to the spin sub-levels of and obtained ^31/^32 = 
1.15 ± 0.05 for one centre and 1.6 ± 0.4 for a second centre. The small 
preference of the decay from the lower energy singlet ^E to the = 
spin sub-level implies that optical spin-polarisation into the rus = spin- 
projection must result from the preferential decay from the ^E rris = ±1 
spin sub- levels to the higher energy singlet ^Ai, rather than the preferential 
decay from ^E to the ^^2 = spin sub-level |181j . This implication can 
be demonstrated more quantitatively by applying the six level rate equation 
model to calculate the ground state spin-polarisation under continuous off- 
resonance excitation |165j 

p = ^ = ^ + ^° ^l-ad (11) 

where P^^ are the populations of the ground state spin sub-levels under con- 
tinuous off-resonance excitation, and a = /C32//C31 and /3 = k^i/k^^ are the 
ratios of the non-radiative decay rates from ^E to the spin sub-levels of ^^2 
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and from the spin sub- levels of to ^Ai, respectively. The approximation 
made in the above expression assumes that 070 ^ a/3 ^ 1. Consequently, in 
order for the spin-polarisation to be ~ 80% (average quoted spin-polarisation 
from the vast range of spin-polarisations reported), a/3 must be ~ 0.2. For 
example, given a ~ 0.8 as determined by Robledo et al, the above indi- 
cates that (3 must be ~ 0.2 in order for Vgs ~ 80%. Thus, it can be seen 
that reasonably high degrees of spin-polarisation can be achieved without 
the extreme differences in the decay rates from the spin sub-levels of "^E to 
^Ai that would be expected if the transitions for one of the spin-projections 
were forbidden. It may be further concluded from the above expression that 
the dependence of spin-polarisation on the ratio of non-radiative transitions 
implies that spin-polarisation will exhibit a significant temperature depen- 
dence. This temperature dependence has been observed by Felton et al [130] . 
who found using conventional EPR techniques that optical spin-polarisation 
decreased from 10-300 K. Spin-polarisation is therefore a reasonably subtle 
effect in NV~, which may depend intricately on temperature, crystal strain 
and other extrinsic factors, thereby explaining the large variance in reported 
values. Hence, spin-polarisation requires a fully developed theoretical model 
for it to be studied in any detail. 

5.3 The spin dynamics of NV 

The spin dynamics of NV~ include the transitions that occur between the 
spin sub-levels of a given triplet state that act to relax any population differ- 
ence between them, as well as the fluctuating interactions that act to dephase 
the spin states. These two distinct processes are characterised (as for any 
spin system) by the spin relaxation rate 1/Ti and the spin homogeneous 
dephasing rate I/T2. A further complication arises when measurements are 
performed on an ensemble of centres, where each centre has different local 
conditions, or when an ensemble of measurements of a single centre are per- 
formed, where each measurement is made in different conditions. For such 
ensemble measurements, the expectation value of a given observable must be 
averaged over the relevant distribution of conditions. This averaging leads to 
inhomogeneous spin dephasing that is characterised by the rate l/T^, such 
that T2* < T2. The spin relaxation and dephasing rates are sensitive to various 
factors including the densities and types of paramagnetic impurities present 
in the vicinity of the NV centre, the applied magnetic, electric and strain 
fields as well as electron-phonon interactions and temperature. Due to the 
many applications employing the spin of NV~, its spin dynamics have been 
extensively studied and a sample of reported values are contained in table fTTl 
Given the scope of this review being restricted to the intrinsic properties of 
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Table 11: Experimental measurements of the NV~ ground state spin relax- 
ation and dephasing times tabulated by reference. The magnetic field data is 
provided in units of Gauss or Tesla if known or operating frequency of EPR 
apparatus if utilised. 



Reference 


T (K) 


Mag. Field 


Ti fms) 


To (us) 


To* (us) 


Redman'^-'' |147i 1161] 
Kpnuedv'^''' 11851 
Kennedy'^-'' |lU8| 
Jelezko^-'^''' [6] 


80 




265 






1 5-100 

X . ±\J\J 


85 GH/ 




32 




300 


35 GHz 




58 




Room 






1-2 




Jelezko^-'' [5] 


Room 






< 6 




Popa^'"'^ |119| 


Room 




100 G 




0.9 
1.2 




GaebeP'" |12U| 


Room 






350 




Childress'-" |122| 
Takahashi^''' |184j 


Room 








1.7 


Room 
40 
20 
1.7 


8 T 
8 T 
8 T 
8 T 


7.7 

3.8x10^ 


6.7 

8.3 

250 




Balasub."'^ [13] 


Room 






1.8x103 




Stanwix^'" |186] 


Room 






600 





"type Ila; ^type lb; isotopically pure type Ila 
'^ensemble; 'single centre; "nanodiamond 



the NV centre, only aspects of the dynamics that provide information about 
such properties will be discussed in detail. 

Van Oort and Glasbeek [183] were the first to measure the dependence of 
the ground state spin relaxation time Ti on an applied magnetic field (refer to 
figure [3T1) . Van Oort and Glasbeek conducted their measurements using NV~ 
ensembles in type lb diamond and showed that the magnetic field dependence 
of the ground state Ti was governed by the resonances of the ground state fine 
structure levels of NV~ sub-ensembles (split by the applied magnetic field) 
and resonances between one sub-ensemble and the electronic spin of the A^^ 
centre. In the regions of these resonances, the NV~ centre can undergo 
energy conserving spin fiip-fiops that redistribute any population difference 
in the spin sub-levels of one sub-ensemble with the other sub-ensembles or 
spins, thereby leading to a relaxation of any initial population difference |148l 
1183] . As demonstrated by figure [HI the ground state Ti is approximately 
magnetic field independent for fields that are away from any resonance. The 
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observations of van Oort and Glasbeek can be extrapolated to single NV~ 
centres by noting that the effects of any resonances with distant NV~ centres 
will be negligible, leaving just resonances with N<j centres to influence Ti. 
The Ns centre is not special in this property as any paramagnetic impurity 
with an electronic g-factor comparable to NV~ will become resonant in some 
magnetic field. This naturally precludes nuclear spins, such as those of ^'^N 
and ^^C, from inducing spin relaxation except in extreme field configurations. 

Away from resonances with paramagnetic impurities, the ground state 
Ti is governed by phonon induced transitions between the spin sub-levels. 
Redman et al |147] were the first to measure the temperature dependence 
(in the range 100-500 K) of the ground state spin relaxation rate 1/Ti of 
an ensemble of NV~ centres. They found that up to room temperature, the 
relaxation rate was well described by the expression 1 /Ti oc co^/ (^e^^/^BT _ 
using w = 63 meV. From the form of this expression and the equivalence 
of the fitted phonon energy with the dominant phonon energy of the opti- 
cal band ~ 67 meV, Redman et al concluded that relaxation occurred via 
a two-phonon process involving strongly interacting modes of energy ~ 63 
meV. Above room temperature, they found that the rate was well described 
by 1/Ti oc T^, which as previously discussed in section |U is indicative of a 
two-phonon Raman process involving the continuum of low energy phonons. 
The Ti measurements of Redman et al were later repeated using a different 
method involving very large axial magnetic fields (~ 8 T) by Takahashi et al 
|184] . In partial agreement with Redman et al, Takahashi et al found that 
1/Ti oc AT + BT^, where A and B are fitting parameters. The linear relax- 
ation process is indicative of direct one-phonon transitions between the spin 
sub-levels, which although present in smaller magnetic fields, are typically 
insignificant |143] . The direct transitions are likely to have gained signifi- 
cant enhancement due to the larger spin transition energies at ~ 8 T, such 
that they became observable in Takahashi et al's measurements. Regardless, 
it is clear that further measurements of the temperature dependence of the 
ground state Ti is required to confirm the phonon relaxation processes that 
are occurring. There is yet to be developed a detailed model of the possi- 
ble phonon relaxation processes that would aid the analysis of any future 
measurements. 

From their first ODMR experiments, van Oort et al [78] have attributed 
the homogeneous dephasing of the NV~ ground state spin to spin-spin in- 
teractions with the electronic and nuclear spin bath formed by other NV~ 
sub-ensembles, the electronic and nuclear spins of the centre and the nu- 
clear spins of the ^^C isotopic impurities. Kennedy et al |108[ 1185] has shown 
that the interactions with the Ng bath is the dominant dephasing mechanism 
in type lb diamond and have correlated the ground state T2 with the N con- 
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Figure 31: The observed axial magnetic field dependence of the NV~ ground 
state (a) relaxation time Ti |183] and (b) homogeneous dephasing rate I/T2 
|121] at room temperature, (c) The observed temperature dependence of 
I/T2 in an axial magnetic field of ~ 8 T |184] . Measurements were made 
in each case using NV~ ensembles in type lb diamond. The axial magnetic 
fields corresponding to the spin level resonances between the ground state 
NV" spin and the N, defect (NV-Pl), two NV" sub-ensembles (NV-NV) and 
the ground state LAC marked in (a). 
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centration present in the diamond. This is supported by the measurements 
of Takahashi et al, who observed that the temperature dependence of I/T2 
was well described by the temperature dependence of the spin flip-flop 
rate, as well as the measurements of Hanson et al |121] . who observed the 
magnetic field dependence of I/T2 to be governed by the resonances with 
in analogy with the magnetic field dependence of 1/Ti (see figure [ST]) . Simi- 
larly, for single centres in type Ila diamond (undetectable N concentration), 
interactions with the ^^C spin bath are the dominant dephasing mechanism 
and T2 has been correlated with the ^^C concentration [T3 |ll05lll87j . Indeed, 
in isotopically pure diamond, the ground state T2 reaches ~1.8 ms [13] and 
the dynamics of the NV~ spin become intricate due to coherent interactions 
with a few proximal ^^C spins, instead of the incoherent interactions with an 
isotropic bath of ^'^C spins present in less pure diamond [188j . Consequently, 
only in isotopically pure diamond does the spin dephasing rate potentially 
have detectable contributions from phonon interactions. No contributions 
have been detected to date. Since the spin dephasing rate is governed in 
most cases by the spin-spin interactions between NV~ and other electronic 
and nuclear spins, which could be equally described by the interactions of a 
canonical triplet spin system with other spins in the lattice, the spin dephas- 
ing rate does not provide any information about the intrinsic properties of 
the centre. 

Apart from macroscopic inhomogeneities in the applied electric, magnetic 
and strain fields and in the preparation of the ensemble of NV~ centres or 
single NV~ centre, the inhomogeneous spin dephasing time Tg* is governed 
by the microscopic inhomogeneities of the local fields present at each NV~ 
centre of the ensemble or the single centre induced by the slow fluctuations 
of other magnetic and electric defects in the lattice. The slow fluctuations 
of the defects alter the local fields at each centre of the ensemble for a given 
measurement or, for the case of repeated measurements of a single centre, 
alter the local fields for each measurement. The inhomogeneous dephasing 
due to paramagnetic impurities and its dependence on the applied magnetic 
field has been well observed and does not in itself yield any information on 
the intrinsic properties of NV~ . However, recently the dependence of the in- 
homogeneous dephasing on the electric and strain fields has been observed by 
Dolde et al [25] in their electrometry demonstration. They observed that Tg* 
increased dramatically when the magnetic and electric fields were such that 
their effects were of comparable magnitudes. Their observation emphasised 
the combined effects of the interactions of the NV~ ground state spin with 
each of the different fields and suggested that the spin's sensitivity to the 
fluctuations of magnetic and electric defects could be controlled by applied 
fields. The combined effects of each of the fields on the NV~ ground state 
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spin is yet to be modeled in any detail and the observations of Dolde et al 
still require a full explanation. 

The relaxation and dephasing of the spin associated with the excited 
state has only been investigated by Fuchs et al [189] in their single centre 
excited state ODMR experiments at room temperature. Within a large axial 
magnetic field (~ 1276 G), Fuchs et al used an intense microwave field to 
manipulate the spin within the lifetime of ^E and implemented Rabi, 
Ramsey and spin-echo pulse sequences. Fuchs et al measured the homo- 
geneous dephasing time of the spin to be T2 ~ 5.8 ns and found that the 
dephasing rate had almost equal contributions from both the decay out of 
^E as well as the spin-conserving phonon transitions between the ^E fine 
structure levels that are responsible for orbital averaging at room temper- 
ature. The spin relaxation time of ^E has not yet been measured. Fuchs 
et al modeled the phonon transitions between the fine structure levels by a 
fluctuating magnetic field, however a rigorous model of the dephasing and 
relaxation processes that also accounts for the orbital averaging and dynamic 
Jahn- Teller effects in ^E is yet to be developed. 

6 Conclusion 

This review has, for the first time, assembled the body of facts that form 
our current understanding of the NV centre. Any future study of the centre 
should seek to resolve the remaining issues concerning the centre, which can 
be summarised as the following: 

1. The microscopic origins of the equilibrium charge state, the photocon- 
version process and spectral diffusion 

2. The temperature variations of the fine structures of the NV~ ground 
and optically excited triplet states 

3. A detailed analysis of the Jahn- Teller effect and its influence on the 
NV~ optical and infrared vibronic bands 

4. The rigorous characterisation of the NV~ optical spin-polarisation and 
readout mechanisms 

It is evident that the resolution of the above issues will alleviate any barri- 
ers in the development of the current applications of NV~ that are related 
to our understanding of the centre and will also enable the proposal and 
development of new applications. 
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The resolution of the first issue will facilitate greater control over the NV 
charge state during fabrication and implementation. The development of a 
model of the common microscopic origins of the charge related phenomena 
requires complementary ab initio and experimental studies of the effects of 
the systematic manipulation of the microscopic distributions of charge donors 
and traps. Possible first experiments could simply involve the observation of 
NV ensembles in samples of varying macroscopic concentrations of isolated 
substitutional nitrogen defects. Alternatively, observations of single ion im- 
planted centres in high purity diamond could be performed as functions of 
the ion implantation parameters. Given the significant crystal volume re- 
quired to realistically model typical concentrations of donors and acceptors, 
the ab initio techniques that have been previously applied to model the NV 
centre are unsuitable. More suitable techniques are likely to be found in 
charge transfer and reaction modeling techniques of chemistry. 

The resolution of the second issue leads to an enhanced understanding 
of how to minimise the decoherence of the NV~ spin due to temperature 
fluctuations in its current applications or how to optimise its future imple- 
mentation as a thermometer. Whilst there exist comprehensive studies of 
the temperature shift of the ground state fine structure, there exists only 
one (ensemble) study of the temperature variation of the excited state fine 
structure from low to room temperatures. Future objectives should thus in- 
clude the more detailed study of the temperature transformation of the NV~ 
excited state fine structure as well as the development of theoretical models 
of the temperature variations based upon the non-radiative dynamics of the 
centre. 

The resolution of the third issue will enable greater control over the op- 
tical transition and its coupling to photonic structures, as well as provide 
insights into how to extend spin-photon entanglement towards room tem- 
perature conditions. Experimental studies should seek to probe the polari- 
sation, stress response and temperature dependence of different components 
of the vibronic bands. These experimental studies should be complemented 
by the systematic application of Jahn- Teller models of increasing complexity 
(i.e. linear vs quadratic vs anharmonic coupling and effective mode vs many 
modes). Other vibronic effects, such the vibronic coupling of energetically 
close electronic states should also be considered during the interpretation of 
the experimental observations. 

The resolution of the last issue will provide the means to enhance the 
preparation and readout fidelity of the NV~ ground state spin for its qubit 
and metrology applications. Given the predominately non-radiative nature 
of the NV~ optical spin-polarisation mechanism, aspects of the mechanism 
can only be indirectly probed via the optical and infrared transitions. Con- 
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sequently, the first steps in characterising the spin-polarisation mechanism is 
to extend the recent works [18H I5U] that focused on measuring the tempera- 
ture dependence of the hfetimes of all electronic levels involved in the optical 
cycle. The inherent spin dependence of the mechanism necessitates the mea- 
surement of the lifetimes of individual fine structure levels. It is therefore 
clear that the second, third and final unresolved issues are intimately linked. 
Theoretical work should thus attempt to identify the correlations between 
the vibronic/ non-radiative aspects of the these unresolved issues in order to 
ensure a consistent resolution of all three. 

As a concluding remark, this review has demonstrated that the NV centre 
is rich in condensed matter physics, including important, but subtle, relativis- 
tic and vibronic effects. Indeed the depth of field of the physics exhibited 
by the NV centre has resulted in the above issues that have remained unre- 
solved for almost 50 years. Yet, having stated this, much of the physics of the 
NV centre can be understood from the simple and intuitive molecular model 
proposed by Loubser and van Wyk [HI |15]. Consequently, the NV centre 
promises to be an excellent pedagogical tool for condensed matter physics 
and model system to aid in the understanding and identification of other 
defects in semiconductors with comparable or better properties for quantum 
technology applications. There is thus great motivation to pursue a more 
complete understanding of the NV centre. 
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Appendix Properties of the C^y point group 

In this appendix, the key properties of the point symmetry group of the NV 
centre are detailed in order to provide a quick reference for those studying 
the centre. The group's symmetry elements, character table, direct-product 
table and Clebsch-Gordan coefficients are provided. The group projection 
operators, irreducible tensor operators and the Wigner-Eckart matrix element 
theorem are defined along with selection rules for the tensor operators. 
The Csy group is a trigonal point symmetry group of order six. The 
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elements of the group are defined witli respect to the trigonal symmetry 
axis and form three classes, being the identity e, the C3 rotations about 
the symmetry axis and the vertical reflections through planes containing 
the symmetry axis. The operations of the six C^^ group elements on the 
structure of the NV centre are depicted in figure [321 Given three classes, the 
group possesses three irreducible representations Ai, A2 and E, being of 
dimensions one, one and two, respectively. 

The character and direct product multiplication tables of the C^v group 
are contained in tables [12] and [131 respectively. Basis elements that trans- 
form as the Ai irreducible representation are invariant under all of the group 
element operations. Basis elements that transform as the A2 irreducible 
representation are invariant under the operations of the group elements of 
the e and C3 classes, but are antisymmetric under the operations of the 
group elements of the ay class. Basis elements that transform as E are non- 
symmetric. Defining the coordinate system such that the z axis is collinear 
with the trigonal symmetry axis, the irreducible representations that the co- 
ordinate components and coordinate rotations transform identified 
in table [12 

Whilst the matrix representations of the one dimensional irreducible rep- 
resentations are as specified by the group character table, there is a degree of 
flexibility in the definition of the matrix representations of the E irreducible 
representation. The E irreducible representation is most simply defined via 
the operations of the group elements on the cartesian components x, y. How- 
ever, the matrix representations of the E irreducible representation may also 
be defined via the operations of the group elements on the spherical compo- 
nents ^± = T^{x ± iy). The matrix representations (including both carte- 
sian and spherical representations of E) of the C^y irreducible representations 
are contained in table [TH 

The Cst, projection operator corresponding to the A row of the j*'' irre- 
ducible representation is defined as |190j 

^S = |Er(^H^?0lA4 (12) 

i 

where Ij is the dimension of the j*'' irreducible representation, Qi is the z*'* 
group element, T^^\gi)x\ is the A*'' diagonal matrix element of the matrix 
representation of the group element in the j'*^ irreducible representation, 
and Pg. is the unitary transformation operator corresponding to the group el- 
ement Qi. Let ip{ be a state that transforms as the A row of the j^^ irreducible 
representation, then satisfies 

-PiM = S.,x6,,i;i (13) 
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Similarly, an irreducible tensor operator Q-'^ that transforms as the A row of 
the j*'* irreducible representation, is an operator that satisfies 

P9MP^' = ^^^'\9^).xQi (14) 

re 

for every group element Qi and every row A. A matrix element {'ip!^\Q'i\4>x) 
of an irreducible tensor operator can be expressed in terms of a set of 
reduced matrix elements {ip^\\Q'^\\(f)'^) and the Clebsch-Gordan coefficients 
( I . . .) of the group via the Wigner-Eckart theorem 



( p 


Q 


r \ 






1 ) 



{i^^'imm (15) 



Importantly, the reduced matrix elements are only dependent on the irre- 
ducible representations of the tensor operator and the basis states and not the 
row of the irreducible representations that they belong to. The set of reduced 
matrix elements of a tensor operator thus represents the minimum set of inde- 
pendent parameters (as determined by symmetry) that completely describe 
the matrix representation of a tensor operator. The Clebsch-Gordan coeffi- 
cients corresponding to the cartesian and spherical representations of the Cs^, 
group are contained in tables [15] and [16], respectively. These Clebsch-Gordan 
coefficients define the cartesian and spherical tensor operator selection rules 
contained in tables [T7] and [TS] respectively. 



Table 12: Character table of the C^y point group |19U] . The irreducible rep- 
resentations that the coordinate components (x, y, z), quadratic coordinate 
components {x^, xy, xz, y"^, yz, z"^) and coordinate rotations {Rx, Ry, Rz) 
transform as are as identified. 
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(d) (Tl (C) (72 (f) ^3 



Figure 32: The operations of the Cst, group elements on the structure of the 
NV centre. The nearest neighbour nitrogen (brown) and carbon (grey) atoms 
to the vacancy (transparent) are as denoted. The operations of the group 
elements result in permutations of the carbon atom labels. The directions 
of rotation (red arrows) about the trigonal symmetry axis (white arrow) of 
the and group elements are depicted in (b) and (c). The planes of 
reflection (red squares) of the ai, a2 and 0-3 group elements are depicted in 
(d), (e) and (f). 



Table 13: The C-^y point group direct product multiplication table [191] . 
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Table 14: The matrix representations of the C^^ group irreducible representations. Both the cartesian (x, y) |19] 
and spherical (,^+, ^J) |192j representations of E are detailed, e = e^^^/^ . 



a. 



0"! 



0"2 
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V3 
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. 
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v/3 
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. 2 








-e* 






Table 15: The non-zero Clebsch-Gordan coefficients corresponding to the 
cartesian matrix representations of the Cs^, irreducible representations |49] . 
The indices j, k denote the element of the matrix defining the Clebsch- 
Gordan coefficient in the basis {x, y}. 



1 1 
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1 1 
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1 1 
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1 3 
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Table 16: The non-zero Clebsch-Gordan coefficients corresponding to the 
spherical matrix representations of the Cs^ irreducible representations. The 
indices j, k denote the element of the matrix defining the Clebsch-Gordan 
coefficient in the basis {^+, ^^}. 
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Table 17: The selection rules for cartesian tensor operators as derived using 
the Clebsch-Gordan coefficients of table [151 The elements of the four tables 
below correspond to the Clebsch-Gordan coefficients of the matrix elements 
('?/'!i|O^|0^) as defined by the Wigner-Eckart theorem. 
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Table 18: The selection rules for spherical tensor operators as derived using 
the Clebsch-Gordan coefficients of table dni The elements of the four tables 
below correspond to the Clebsch-Gordan coefficients of the matrix elements 
(■01^10^10^) as defined by the Wigner-Eckart theorem. 
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